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' We review the model building aspects for keV sterile neutrinos as Dark Matter candi- 

dates. After giving a brief discussion of some cosmological and astrophysical aspects, we 
first discuss the currently known neutrino data and observables. We then explain the 
purpose and goal of neutrino model building, and review some generic methods used. 
Afterwards certain aspects specific for keV neutrino model building arc discussed, be- 
fore reviewing the bulk of models in the literature. We try to keep the discussion on a 
pedagogical level, while nevertheless pointing out some finer details where necessary and 
useful. Ideally, this review should enable a grad student or an interested colleague from 
cosmology or astrophysics with some prior experience to start working on the field. 
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1. Introduction 

The Dark Matter problem has been a long standing one in physics: since Oort 
in 1932 (by studying stellar motions 1 ) and Zwicky in 1933 (by studying galaxy 
clusters 2 ) independently concluded its existence, Dark Matter (DM) has puzzled 
the scientific community. Even though we know that it must exist, 3 we do not 
know much about its true Nature. While alternative explanations such a Massive 
Compact Halo Objects (MACHOs) 4 or Modified Newtonian Dynamics (MOND) 5 
now seem to be strongly disfavored by the observation of the Bullet Cluster, 6 we 
tend to believe more and more that DM consists of some new, yet unknown, type 
of elementary particle. 7 Depending on its velocity, one classifies particulate DM 
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into Cold (CDM), Warm (WDM), a or Hot (HDM) Dark Matter. While HDM is 
constrained by cosmological structure formation to make up at most about 1% of 
all DM, 9, 10 it is not clear which of the other two possibilities is the correct choice. 
Even mixed DM scenarios are possible. 11, 12 

An equally puzzling topic in elementary particle physics is the Nature of the 
neutrino: after its postulate by Pauli 13 in 1930 and its first detection by Cowan and 
Reines 14 in 1956, we have gone a long way to measuring all light neutrino mixing 
angles 15-17 and narrowing down its mass range to be below the eV-scale. 3 ' 18-22 
Nevertheless, not only have we not been able to experimentally pin down the actual 
neutrino mass, its apparent smallness and also the oddly large mixing angles look 
as if some unknown pattern was behind them. While at the moment we have no 
fully accepted explanation for this pattern, we have yet been successful in building 
models that can explain the measured patterns and at the same time lead to testable 
predictions. 

The present review deals with a crossover topic that relates both fields described 
above: model building for keV sterile neutrino Dark Matter. In this approach, typi- 
cally the first-generation sterile neutrino Ni is the lightest and plays the role of the 
DM particle. The minimal setting yielding such a pattern is the so-called neutrino- 
minimal Standard Model (j/MSM), 23 a framework where the SM is simply extended 
by three right-handed neutrinos with a very specific mass pattern. This framework 
can accommodate for a vast variety of phenomena, such as neutrino oscillations, 24 
the baryon asymmetry of the Universe, 25 ' 26 or the DM problem. 23 However, the 
i/MSM does not contain any explanation for the required mass pattern, i.e., it has 
no mechanism enforcing it. This is the reason why more specific models are required, 
which contain a suitable mechanism that can to some extent explain the appearance 
of the required mass pattern. We will review the current models on the market, and 
explain their respective advantages and disadvantages whenever possible. The key 
point is that a sterile neutrino with a mass of a few keV could be a very suitable 
(typically warm) DM candidate. If this is the case, we will have a direct relation 
between Dark Matter and neutrinos. Connected to that, limits from either sector 
will constrain the other. 27 This raises the question of whether it is actually possible 
to find complete models that can at the same time lead to a keV sterile neutrino 
and to successful active neutrino mixing. It will turn out that this is indeed possi- 
ble. However, due to constraints from the DM sector one has to take into account 
different and new aspects compared to ordinary model building for light neutrinos. 
Generically the models reviewed here are quite strongly constrained, since they have 
to explain data from various sectors. On the other hand, this is exactly what makes 
keV sterile neutrino DM so interesting: if such a neutrino remains to be a good 
DM candidate, it will be necessary that astrophysicists, cosmologists, and particle 
physicists collaborate on the topic in the future, in order to prove or disprove the 
role of the keV sterile neutrino as DM. 



a One might also find references to Cool (CoolDM) Dark Matter, see e.g. Rcf. 8. 



26,2013 3:10 WSPC/INSTRUCTION FILE keV'Merle 



4 Alexander Merle 

While even the astrophysical and cosmological aspects of keV sterile neutrino 
DM had hardly been discussed for much more than one decade, the field of model 
building for keV sterile neutrinos is even younger and less developed, which suggests 
the question why at all one could need a review on the subject at this stage. Fur- 
thermore, the current review is written by a person whose main fields of research 
are neutrino and DM phenomenology and model building which will, naturally, lead 
to a somehow biased view. However, since the field strongly requires experts from 
different disciplines to collaborate in order to investigate it thoroughly in all its 
aspects, it is certainly useful to have a review at hand which explains the goals 
and methods of one of the related disciplines in a preferably pedagogical matter. 
While very good overviews of the astrophysical and cosmological aspects of keV 
sterile neutrinos already exist, 28 ' 29 the present manuscript comprises an attempt to 
achieve a similar work for the model building aspects of keV neutrinos. Ideally, the 
text should serve as a travel guide for astrophysicists and cosmologists who have 
been puzzled for long by the nitty-gritty pieces of particle physics model building. 
At the same time, while trying not to be prejudiced, the review is supposed to 
provide the reader with the argumentations used by model builders when consid- 
ering certain aspects of a model as "good" or "bad" , or rather "advantageous" or 
"disadvantageous" . On the other hand, particle physics graduate students should be 
equally able to profit from the text which gives a pedagogical overview of the matter 
without going into formal details whenever they can be avoided. Finally, even one 
or the other more experienced researcher in particle theory might be interested in 
a concise summary of the topic or in a motivation for why to study keV neutrinos. 
The sincere hopes of the author are that these multiple needs are covered by the 
present review. 

The text is structured as follows. After shortly discussing some astrophysical and 
experimental aspects of keV sterile neutrinos in Sec. 2, we review the experimental 
status and the construction of models for active neutrinos in Sec. 3. We then give in 
Sec. 4 a brief discussion of some aspects which are specific for keV sterile neutrinos 
and do not usually play a role for light neutrinos. The central chapter of the review 
is Sec. 5, where most of the models currently present in the literature are discussed 
in various detail, depending on feasibility of and suitability for a pedagogical dis- 
cussion. We finally conclude in Sec. 6. Some technical details on the neutrino mass 
matrix are discussed in Appendix A. 

Before starting the discussion, it is worth to point out that the reader should 
at least be familiar with the basics of the fields. Some knowledge on the Standard 
Model (SM) of Elementary Particle Physics and on the aspects of quantum field 
theory required to understand the SM is absolutely necessary to understand the 
text. Furthermore, some very basic group theory (at a level that even the author 
can understand) is used without sufficient explanation, while the more complicated 
aspects are explained in more detail. Finally, since the text focuses on the particle 
physics aspects of keV sterile neutrinos, an elementary knowledge of cosmology is 
not strictly necessary but strongly recommended. To finish with a comment on the 
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notation, throughout this review we will make use of natural units: H = c = 1. 
2. Astrophysical and other Aspects 

While the question about the nature of DM is mainly a particle physics problem, 
many of the properties of DM are known from astrophysics. Although we cannot 
review all the astrophysical requirements here, we want to at least mention them and 
give a flavor of their importance. We also give some references for further reading, 
regarding additional implications of or probes for keV sterile neutrinos. 

2.1. Warm or cold?!? 

HDM is excluded as the dominant DM component by cosmological structure for- 
mation, 9,10 since the corresponding top-down formation (i.e., large structures form 
first) leads to contradictions with observations. However, the situation of WDM 
and CDM is not quite as clear. While historically mostly CDM has been consid- 
ered, due to the variety of possible candidates in theories with, e.g., supersymmetry 
or extra dimensions, at the moment we unfortunately do not seem to have any 
solid direct evidence for such a DM candidate. Even worse, recent experiments such 
as XENON 30 ' 31 cut more and more chunks from the allowed parameter space, and 
some CDM candidates may at some point start being in trouble. 

WDM, on the other hand, has also been investigated for some time by the 
structure formation community: while pioneering studies appeared already more 
than one decade ago, 32 ' further collaborations have formed by now which inves- 
tigate the topic. 12,34 One particular focus of the simulations is the structure at 
relatively small scales, 35-39 since this is the main point distinguishing WDM from 
CDM. Thus smaller galaxies, and in particular the so-called dwarf satellite galaxies, 
are considered to be the key to distinguish the two types of DM. Indeed, WDM sim- 
ulations tend to lead to a smaller number of dwarf satellites, 34 which seems to be 
in coincidence with the fact that we seem to observe considerably less dwarf satel- 
lites in space than predicted by CDM simulation. This is often called the missing 
satellite problem 40 ^ 2 in the literature. However, there could also be astrophysical 
reasons for this, e.g., it may be that supernova explosions occur in the formation 
of dwarf satellites at a high enough rate to simply blow away much of the visible 
material, 43 ' 44 such that the dwarf satellite galaxies appear much fainter than larger 
galaxies. In any case, there is agreement that the core profiles of dwarf satellites 
are the key to distinguish between CDM and WDM, and hence the dwarf satellites 
which have already been observed should be investigated in more detail (see, e.g., 
Refs. 34, 45, 46, 47, 38). Note also that, on the more observational side, model- 
independent surveys 48 and data analyses 49 seem to point towards a keV-mass DM 
particle. 

To give a flavor of the lower bounds one can obtain on the keV neutrino mass, 
the most model-independent limit one can derive originates from the Tremainc- 
Gunn bound. 50 The principle idea is that the averaged phase base density of a 
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fermionic DM candidates cannot be smaller that of a degenerate Fermi gas. When 
applied to the case of keV sterile neutrinos, the typical lower bound one obtained 
on the DM mass is around M > 1 keV. 51,52 Stronger limits can be obtained when 
the production mechanism is taken into account, which results in different lower 
bounds ranging from 1.6 to 10 keV, see e.g. Refs. 53, 54, 55, 56. 

Let us end this section by stressing that, although often used more or less equiv- 
alently in the literature, a keV mass of the DM particle does not necessarily mean 
that it is WDM. The decisive point is the velocity profile rather than the mere mass 
(e.g., axions have tiny masses but are typically CDM 57 ~ 60 ). Nevertheless, it is true 
in many cases that keV mass particles, and in particular keV sterile neutrinos, in- 
deed turn out to be WDM. But the exact velocity profile depends on the production 
mechanism under consideration, and sometimes on subsequent effects modifying the 
expansion rate and hence the cooling of (parts of) the Universe. 

2.2. Production mechanisms 

A necessary condition for any DM candidate particle is to participate in a mecha- 
nism which can produce enough of this particle to make up all or at least a significant 
part of the amount of the DM in the Universe. This is non-trivial, since the require- 
ment of having not too much and not too little DM can often only be fulfilled in 
a very narrow parameter range. The most generic DM production mechanism is 
the so-called thermal freeze-out, which is discussed in many textbooks (see, e.g., 
Ref. 61). However, for this production mechanism to be effective a particle needs an 
interaction strength comparable to that of weak interactions, which is gencrically 
difficult for a sterile neutrino. On the other hand, if the SM gauge group is extended, 
thermal freeze-out might be revived, but this will lead to other difficulties. An alter- 
native thermal production mechanism is freeze-in production, 62 which is well-suited 
for extremely weakly interacting particles. Finally, non-thermal production is also 
a possibility. 

We will now briefly discuss the most generic mechanisms used to produce keV 
sterile neutrino DM. 

2.2.1. The Dodelson-Widrow mechanism 

The easiest and most natural possibility is the so-called Dodelson-Widrow (DW) 
mechanism. 63 The idea, based on Refs. 64, 65, 66, 67, 68, 69, 70, 71, is that al- 
though the interaction strength of keV sterile neutrinos is too small for them to be 
in thermal equilibrium with the plasma, they could nevertheless be produced by 
the plasma. This is because, as we will see in Sec. 3.2.1, sterile neutrinos are typ- 
ically not absolutely sterile but instead have tiny admixtures to active neutrinos. 
Although these admixtures are not large enough to keep the keV sterile neutrinos 
in thermal equilibrium with the plasma, from time to time they are produced in 
the plasma by processes which mainly generate active neutrinos, but which can also 
produce stcrilcs by their admixtures. Thus, a certain amount of keV sterile neutrinos 
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gradually builds up, and the corresponding interaction strengths and densities arc 
small enough so that the keV neutrinos do not annihilate again. This is very similar 
to what happens in freeze-in. 62 While the details of this mechanism can be quite 
tricky (e.g. due to hadronic contributions to the sterile neutrino production 72 ), the 
general idea is nevertheless very easy to understand. 

Despite its simplicity, this mechanism is by now known to be excluded by the 
strong X-ray bound (cf. Sees. 2.3 and 4.1) in the minimal version, i.e., if no primor- 
dial lepton asymmetry, which must arise from further new physics, is present. 29, 73 
Hence, the DW mechanism alone is not enough to produce a sufficient amount of 
DM. On the other hand, the DW contribution to the DM relic density is unavoidable 
as long as the active-sterile mixing is not completely switched off. Thus, many other 
production mechanisms should actually be accompanied by an additional amount 
of keV sterile neutrino production by the DW mechanism. 

2.2.2. The Shi-Fuller mechanism 

An alternative possibility is the so-called Shi-Fuller (SF) mechanism, 8 which can 
arise in addition to the DW mechanism. The basic idea is that a significant primor- 
dial lepton- antilcpton asymmetry in the early Universe can act onto active-sterile 
neutrino transitions in a way similar to a background made of ordinary matter acts 
onto ordinary neutrino oscillations. In such an environment it is well known that 
resonant flavor transitions can occur in the presence of a certain matter density, 
which is known as the Mikheev-Smirnov-Wolfenstein (MSW) effect. 74-78 Such tran- 
sitions can, in the early Universe, produce a large amount of keV sterile neutrinos at 
a very specific temperature, in case that a suitable lepton asymmetry is present. 79 
While this mechanism cannot stand alone, as it will be necessarily accompanied 
by a DW contribution, its great benefit lies in the fact that the neutrino spectrum 
is changed: the non-thermal component produced by the resonant transition leads 
to an overall cooler DM spectrum than for pure DW production. 8 This makes the 
keV sterile neutrinos in the Universe more similar to ordinary CDM, which helps 
to avoid problems associated with a too warm spectrum. 

2.2.3. Scalar decays 

Another option for producing keV sterile neutrinos non-thermally is by the decay of 
scalar singlets 28, 56,80 which could, e.g., freeze-out by Higgs portal interactions and 
decay at a temperature similar to their mass (or, rather, at their freeze-out temper- 
ature which is in turn similar to their mass). One case that has been particularly 
studied is the one of this scalar particle being the infiaton. 81-83 This possibility is 
particularly attractive from the point of view of solving the additional problem of 
inflation, while a general singlet scalar is hardly constrained. Moreover, this mech- 
anism also leads to a somewhat cooler spectrum, 56 which can help to evade certain 
bounds. 
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In addition it is worth to stress the versatility of scenarios with scalar produc- 
tion: not only do they have an immediate connection to the scalar sector of the SM 
and theories beyond, leading to interesting collider signatures, 84 a non-zero vacuum 
expectation value of the scalar field can also immediately lead to Majorana neu- 
trino masses 82 and to other lcpton number violation phenomenology. 85 From the 
particle physics point of view, this production mechanism might be one of the most 
interesting ones. 

2.2.4. Thermal overproduction with subsequent entropy dilution 

Finally, it could also be that the kcV sterile neutrinos are only sterile with respect to 
SM interactions, but are non-trivially charged under the actual (larger) gauge group 
which breaks down to the SM at a low enough temperature. The prime example 
of such a situation is Left-Right symmetry, 86 - 90 SU(S) C x SU(2) L x SU(2) R x 
[/(1)b_l, where all right-handed fcrmions including neutrinos are singlets under 
SU(2)l, but charged non-trivially under SU(2)r. Hence, at temperatures above 
the breaking SU(3) C x SU(2) L x SU(2) R x U{1) B -l -> SU(S) C x SU{2) L x U(1) Y , 
the keV "sterile" neutrino can actually be in thermal equilibrium with the thermal 
plasma in the early Universe, and at some point it can freeze-out like nearly all the 
other particles. 

This idea has been applied to the case of keV sterile neutrinos in Refs. 53, 91, 92, 
93. However, it turns out that, just as for the freeze-out of SM neutrinos if they had 
masses of a few keV, the keV steriles would freeze-out at a temperature much larger 
than their mass, and hence their abundance would not be suppressed. Accordingly, 
they would overdose the Universe by comprising far too much DM. While this is 
seemingly a big problem, it can be compensated to some extent by a subsequent 
dilution by entropy production 94 from the decay of a frozen-out non-relativistic 
species that temporarily dominates the energy density of the Universe. Typically, 
the second and third generation sterile neutrinos N^^ do this job, as their decays are 
generically much faster than those of the Ni due to their larger mass, cf. Sec. 2.3. 
Indeed, such a mechanism can lead to consistency with most bounds, 53,91 but one 
has to mention that such settings are nevertheless under a certain amount of tension 
when all the data is put together. One particularly hard bound is the lower value 
of the "reheating" temperature before the event of Big Bang Nucleosynthesis, 95 ' 96 
which can drastically reduce the allowed parameter space even for relatively general 
settings. 97 

2.3. The X-ray bound 

Probably the most important bound in practice on models with keV sterile neu- 
trinos is their radiative decay, N% — > ^. 98 > 99 This decay proceeds via a 1-loop 
diagram, which looks very similar to typical 1-loop diagrams for the lepton flavor 
violating decay \i — \ ey. 10 °- 103 As generic for such a decay, if the final state lepton 
is practically massless, the decay rate is proportinal to the fifth power of the initial 
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state mass. 104 Accordingly, a decay Ni — > is proportional to Mf , where M, is the 
mass of Ni. If out of three sterile neutrinos the first generation fermion Ni is the 
keV neutrino, while N2,3 are considerably heavier, this implies that the two heavy 
neutrinos decay extremely quickly, whereas the keV neutrino is practically stable 
even on cosmological scales. Otherwise it could not play the role of DM. 

Nevertheless, if enough of the keV neutrinos are present in the Universe, some of 
them will decay and produce a certain amount of photons. Since the light neutrino 
mass is practically zero compared to the mass of the keV neutrino, the photon 
resulting from the decay is practically monoenergetic. Accordingly, the smoking 
gun signature to search for is a monoenergetic X-ray line which can be searched for 
by satellite experiments, see e.g. Refs. 105, 106, 107, 108, 109, 110, 111, 112, 113, 
114, 115, 116, 117, 118. 

Turning the argumentation round, the non-observation of this X-ray line pro- 
duces strong bounds on the active-sterile mixing and hence on the full neutrino 
mass matrix. However, we will postpone the discussion of this constraint to the 
dedicated Sec. 4.1. The reason is that we will first need to introduce some more 
aspects of neutrino physics before we can appreciate the effect of this bound on a 
concrete model, which will be done in the next section. 



2.4. Other aspects of keV sterile neutrinos 

There are many further interesting aspects of keV sterile neutrinos, which are too 
numerous to all be mentioned here. However, to give at least an idea of the many 
existing interrelations, we mention a few. 

First, there are astrophysical implications different from structure formation or 
DM production, where keV sterile neutrinos could be interesting: for example, they 
could help to explain the observed motions of pulsars ("pulsar kicks"), 7 ' 56,119-121 
they could influence star formation, 122 or their decays could lead to the production 
of X-rays causing a partial reionization of the Universe. 123 

One could also think of detection prospects in ground-based experiments: while 
it is unlikely to experimentally probe their tiny interactions with matter, 124 the 
indirect detection of the lightest sterile neutrino N\ in /3-transitions 125-129 or of 
their heavier brothers N23 in pp-collisions and meson decays 130 is at least not 
hopeless. 



3. Neutrino Model Building 

In this section, we will review the model building aspects of light, i.e., eV-scale active 
neutrinos. Although the main focus of this review is on keV-scale sterile neutrinos, 
reproducing the correct low energy neutrino data is a necessary requirement for any 
model to be built. Furthermore, the basic principles of model building can be nicely 
illustrated for light neutrinos. 
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3.1. Neutrino masses and mixings 

After having been theoretically predicted by Pauli 13 in 1930 to rescue the conserva- 
tion laws in nuclear beta decay, the neutrino was detected for the first time in the 
famous Cowan-Reines experiment. 14 At that time, the neutrino has been assumed 
to be essentially massless. On the other hand, a small mass would not have had a big 
effect in any experiment, and could hence not be excluded. However, as proposed 
by Pontecorvo first for meson-oscillations 131 and later for neutrino-antineutrino sys- 
tems, 132 even a tiny neutrino mass can still have an observable effect if the neutrino 
mass basis is not equal to the flavor basis. The resulting phenomenon is known as 
neutrino oscillations, and it essentially means that a neutrino produced in one defi- 
nite flavor, say v e , can change its identity while traveling from one place to another, 
to be finally detected as, e.g., v^: neutrinos are mixed. An excellent although slightly 
outdated review on this subject can be found in Ref. 133. Note that, since neutrino 
oscillation probabilities only depend on mass square differences Am^ = mf — m 2 , 
neutrino oscillations cannot give any information about the actual scale of neutrino 
masses. However, the observation of at least one non-zero Ato 2 ^, i.e. of neutrino 
oscillations, immediately implies that not all neutrinos can be massless. 

3.1.1. Neutrino mixing and neutrino oscillations 

What neutrino mixing means in practice is that, e.g., an electron neutrino v e (i.e., 
the neutrino which couples to the electron via VT-bosons) does not have a definite 
"electron neutrino mass" , but it is rather a superposition of (usually three) mass 
eigenstates Vi, each of which has a well-defined and definite mass m.^. This implies 
that a neutrino produced with a definite flavor will change its identity when prop- 
agating, since the different masses of the different mass eigenstate components will 
lead to different quantum mechanical phases. In the two-flavor approximation, i.e., 
when having two weak interaction eigenstates (14, i^) and two mass eigenstates 
V2), the corresponding oscillation formula for a neutrino produced as a pure v e 
is given by: 133 



where E is the neutrino energy, L is the distance between production and detection 
of the neutrino, and Am 2 = m 2 — m 2 . Indeed, if m\ = mi (and hence in particular if 
m l,2 = 0), the transition probability from Eq. (1) is zero. Furthermore, if the mixing 
angle 9 which describes the mismatch between the two bases is zero, then we also 
have P = 0. In that case, the flavor would completely determine the mass and any 
quantum mechanical uncertainty would vanish, thereby killing the interference-like 
oscillation. More information on the correct treatment of neutrino oscillations in 
the context of quantum mechanics and quantum field theory can, e.g., be found in 
Refs. 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145. 

Formally, 3-flavor neutrino mixing can be described by a rotation of the ba- 




(1) 
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sis, and the corresponding mixing matrix 146 is known as the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) matrix U = C/pmns : 

or v a = U a iVi for short, in a basis where the charged lepton mass matrix is diag- 
onal. Explicitly, the PMNS matrix can be written in terms of three mixing angles 
(012, 013, #23), one Dirac CP-violating phase S, and two Majorana CP-violating 
phases (a, /3), the latter two having trivial values or 7r in the case of Dirac neu- 
trinos: 

/ C12C13 S12C13 si3e~ l5 \ 

U = -S12C23 - ci 2 S23Si3e l<5 C12C23 - si 2 S23Si 3 e iA S23C13 diag(l,e iQ ,e l(/3+,5) ), 

\ S12S23 - Ci2C 2 3Sl 3 e 1 ' 5 -C12S23 - Sl2C 2 3Sl3e i<5 C23C13 / 

(3) 

where Cy = cos0jj and = sin0jj. If these mixing angles were small, we could 
approximate the flavor by the mass eigenstates, similar to the quark sector: there, 
a weak interaction eigenstate down-quark d! is not very different from the corre- 
sponding strong interaction (and mass) eigenstate d, a fact that is reflected in the 
quark mixing angles being relatively small. 147 However, experimentally - and in- 
deed surprisingly at first sight - the neutrino mixing angles have been measured to 
partially be very large: 15 for a normal mass ordering, mi < 7772 < 1713, the mixing 
angles are given by 

sin 2 6»i2 = 0.320 (0.303-0.336, 0.27-0.37), 

sin 2 6»i3 = 0.0246 (0.0218-0.0275, 0.017-0.033), 

sin 2 6*23 = 0.6 1 3 (0.400-0.461 & 0.573-0.635, 0.36-0.68), (4) 

while for inverted ordering, 777,3 < m i < m 2, we have 

sin 2 0i2 = 0.320 (0.303-0.336, 0.27-0.37), 

sin 2 0i3 = 0.0250 (0.0223-0.0276, 0.017-0.033), 

sin 2 023 = 0.600 (0.569-0.626 , 0.37-0.67), (5) 

where we have quoted the best-fit value (la region, 3cr region). Note that the global 
fit data quoted in Eqs. (4) and (5) includes the very recent (2012) measurements of 
the previously unknown mixing angle i3 by the Daya Bay, 148 RENO, 149 and Dou- 
ble Chooz 150, 151 collaborations. Other mass ordcrings different from normal and in- 
verted ordering are not possible, since the MSW effect 74-78 for solar neutrinos, 152-158 
supported by reactor neutrino data, 159 ' 160 forces the mass square difference Am^ = 
777^ - 777 2 to be positive, i.e., Am|i = 7.62 (7.43-7.81, 7.12-8.20) • 10~ 5 cV 2 . 15 

Finally, measurements of atmospheric neutrinos 161 together with experiments 
detecting neutrinos that are artificially produced in accelerators 162-166 can be used 
to obtain |Ato|i| = 2.55 (2.46-2.61, 2.31-2.74) • 10" 3 eV 2 (normal ordering) and 
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(Am^l = 2.43 (2.37-2.50, 2.21-2.64) • 10~ 3 eV 2 (inverted ordering). 15 At the mo- 
ment, we do not know the sign of Ara^, i.e., if mj < m 3 or Hi\ > m 3 , which means 
that the mass ordering has not yet been determined. 

Note that we could, of course, have used any of the other currently available 
global fits to neutrino mixing data, 16,17 but for dcfinitcness we had to decide for 
one and have taken Ref. 15. 



3.1.2. Limits on neutrino masses 

Even though we have considerable knowledge about the two mass square differences 
Atrial and lAmjjjJ, this is not enough to draw conclusions about the absolute neu- 
trino mass scale. We do, however, have experimental bounds on certain observables. 
At the moment, we know more or less three realistic probes of the absolute neutrino 
mass scale, which are its kinematical determination in laboratory experiments on 
single beta decay, the measurement of the effective mass in neutrino-less double 
beta decay, and the cosmological determination of the sum of all neutrino masses. 
Although several present and future neutrino detectors might also be interesting 
in the context of supernova neutrinos (see, e.g., Rcfs. 167, 168), they necessarily 
would need a supernova to happen in our Galaxy in the first place, which is not 
guaranteed. Alternative proposals to get information on the absolute neutrino mass 
(e.g., Rcfs. 169, 170, 171, 172, 173) seem even less realistic. 

In single beta decays, one tries to use nuclear reactions such as (Z, A) (Z + 
1, A) + e~ + 9 e to determine the neutrino mass, which is the simplest and most 
model-independent way to get information on the neutrino mass. The trick is to use 
only the small fraction of all decays in which the electron carries away nearly all the 
kinetic energy, i.e., the emitted neutrino is essentially at rest. Since current and near 
future experiments are not sensitive enough to resolve the three mass eigenvalues 
separately, the corresponding observable arising from an incoherent sum over all 
possible final states, and it is called the effective electron neutrino mass Hip, defined 
by 174 

2 _ 222, 222, 22 r n\ 

m p = m lCl2 c 13 + m 2 s 12 c 13 + m 3 s 13 . (6) 

Note that this mass is really just an effective observable and in no sense related 
to something like a mass of the v e , since that cannot even exist as argued above. 
Although such single beta decay experiments require an enormous amount of pre- 
cision, there are impressive limits by past experiments using tritium: the MAINZ 
experiment measured m 2 p = (-1.2 ± 2.2 sta t. ± 2.1 sys .) eV 2 , 18 and TROITSK mea- 
sured m 2 = (—1.0 ± 3.0 s tat. ± 2.1 sys .) eV 2 , 19 both amounting to an upper limit 
of about 2 eV on the absolute neutrino mass scale. The next generation experi- 
ment KATRIN 175 will have a sensitivity on mp that is better by about one order 
of magnitude. An alternative approach will be taken by the MARE experiment, 176 
using rhenium, which has a lower Q value and hence an improved sensitivity on 
the neutrino mass, but one pays the price of the transition of interest being a first 
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forbidden decay. 

One of the biggest questions in neutrino physics is about their actual nature - 
are the neutrinos Dirac particles, i.e. distinct from their antiparticles, v c ^ e^v, 
or are they Majorana particles, i.e. identical to their charge conjugates up to a 
phase, v c = e^vl If the latter was the case, then an antincutrino v e emitted in a 
f}~ decay (Z, A) — > (Z+ 1, A) + e~ + D e could be absorbed as neutrino in an inverse 
/? decay v e + (Z + 1, A) — > (Z + 2, A) + e~ . If these two transitions happen in one 
single nucleus, with the neutrino being purely virtual, the resulting net reaction is 
called neutrino-less double beta decay (0vf3/3), (Z,A) — > (Z + 2, A) + e~ + e~ . There 
are many more interesting connections of Oz//3/3 to particle physics, as can be seen 
from the very detailed and informative review by Rodejohann, Ref. 177. A review 
that is dedicated a bit more to the nuclear physics side can be found in Ref. 178. 
Recent overviews of the on-going and future experimental activities are provided 
by Refs. 179, 180. 

It is not a priori clear if light neutrino exchange is indeed the dominant contri- 
bution to 0^/3/3 (see, e.g., Refs. 181, 182, 183, 184, 185, 186, 187, 188, 189, 190 for 
discussions of situations where this is not the case). However, if it does dominate, 
then the decay rate is proportional to the square of a quantity called the effective 
mass \m ee \ given by 

\m ee \ = \m lC 2 12 c 2 13 + m 2 s 2 12 c 2 13 e ia + m 3 s 2 13 e il3 \. (7) 

Note that a and /3 are the Majorana phases from Eq. (3), which did not appear 
in the observable m|. This is a reflection of lepton number being violated in the 
reaction (Z,A) — > (Z + 2, A) + e~ + e~ . Note that, contrary to some statements 
in the literature, the effective mass \m ee \ does not depend on the Dirac phase S, 
since this phase can be absorbed in a redefinition of the neutrino mass eigenstate 
v 3 , cf. Refs. 191, 192. In contrast to m|, \m ee \ cannot be split into different parts 
even with an infinitely accurate experimental setup, since the contributions of the 
different neutrino mass eigenstates are purely virtual, and hence the total amplitude 
is given by a coherent sum over partial amplitudes. 

In order to be experimentally observable, one needs a nucleus in which single 
beta decay is kinematically forbidden, as otherwise this process would completely 
dominate the decay. Even then the process of double beta decay with the emission of 
two neutrinos, (Z,A) — > (Z + 2, A) + e~ + e~ + D e + i> e , cannot be forbidden, and one 
has to search for the smoking gun signature of the sum of the two electron energies 
being equal to the Q-value of the decay. Among the isotopes investigated so far and 
in the future are 76 Ge (GERDA, 193 Heidelberg-Moscow, 194 IGEX, 195 Majorana 196 ), 
130 Te (CUORE, 197 CUORICINO 22 ), or 136 Xe (EXO 21 ). Alternative decay modes 
are mainly investigated in COBRA. 198 Up to now, 0i//?/3 has not been observed 13 , 
and typical lower limits on the half-lives are of the order of 10 23 to 10 25 years, which 



b The claim made in Ref. 199 is under strong pressure, if not even excluded, by Ref. 21. 
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translates, modulo nuclear physics complications, into upper limits on \m ee \ of a few 
tenths of an cV, cf. Rcf. 177. 

Finally, note that a positive signal of 0vf3j3 would unambiguously show that lep- 
ton number is violated and that neutrinos arc hence Majorana particles. 200 However, 
the corresponding mass contribution generated by the resulting Butterfly diagram 
would be at most about 10 -24 eV or in some cases even vanishing at 4-loop level, 201 
and hence by far too tiny to explain the minimum neutrino mass scale enforced by 
the measured value of |Am|i|. 

A third way to get information on the absolute neutrino mass scale is by cos- 
mological observations, where neutrinos also play an important role. 202 Using basic 
cosmological arguments, one can determine the ratio of the neutrino temperature in 
the Universe after e + e~ annihilations to the "reheated" photon temperature. This 
allows to relate the sum E of all light neutrino masses to their relic density Vl v h 2 : 



if one now uses the observations of the cosmic microwave background (CMB) to 
determine the total matter density f2 ma t, as well as the baryon density and the 
dark matter relic density f2DM c , one can determine the neutrino relic density as f2„ = 
fi ma t — ilb — £!dm- 203 The newest analysis of the Wilkinson Microwave Anisotropy 
Probe (WMAP) 9-year data then yields a limit of S < 1.3 eV @ 95% C.L., 204 which 
can be improved to 0.44 eV when including data from other CMB-measurements, 
from baryon acoustic oscillations, and from supernova redshift surveys. Apparently, 
these limits depend strongly on which data set and which type of statistics are used. 
Note that there has been a recent hint for a non-zero value E = (0.32 ± 0.11) eV 
by the South Pole Telescope (SPT) collaboration. 205 However, the similar Atacama 
Cosmology Telescope (ACT) 206 currently does not seem to confirm this observation. 
Accordingly, we will have to wait for new data in order to confirm or refute this 
hint. 

The aim would be to combine the cosmological value of E with laboratory data, 
which could positively influence each other. 207 Ideally, we would therefore like E to 
be the sum of only the three active neutrino masses, 



However, E actually denotes the sum over all light neutrino species, which could 
also include sterile neutrinos. 208 In case that one erroneously assumed Eq. (9), but 
the measured parameter E is in fact dominated by other contributions {e.g., by 
systematic errors we might not even be aware of), a naive analysis could lead to 
wrong limits, or even to a fake "measurement" of the neutrino mass scale. 209 

Summing up, we have, in particular within the last decades, gained considerable 
knowledge on the neutrino masses and mixing parameters. Yet we are puzzled by 



E = 94 cV ■ n„h 2 . 



(8) 



E = mi + m 2 + m.3. 



(9) 



c To be precise, the standard cosmological analysis actually uses the ACDM model, i.e., the cos- 
mological standard model including Dark Energy and Cold Dark Matter. 
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their values, and in particular by the smallness of their masses and by the large 
mixings compared to the quark sector. To the best of our knowledge, this seems to 
suggest a structure, or symmetry principle, behind these parameters. It is exactly 
at this point where neutrino model building comes into the game. 

3.2. Model building for light neutrinos 

The goal of model building in neutrino physics can be summarized in one sentence: 
we aim to find a particle physics explanation for the sizes and values of neutrino 
masses and leptonic mixing angles. 

Not every approach is powerful enough to explain both these aspects, so it 
makes sense to divide the models on the market into mass models, which give an 
explanation for the tininess of the active neutrino masses, and flavor models, which 
can explain the mixing pattern in the leptonic sector. d Typically, mass models do 
not predict an exact value for the mass, which is intrinsically difficult in a quantum 
field theory, but rather some hierarchy among certain masses, e.g., why the light 
neutrinos should be much lighter than the other SM fermions. Flavor models, in 
turn, often predict very specific values for certain mixing angles, however at the 
prize of introducing practically unobservable high energy sectors. On top of that, 
the magnitude of mixing angles is intrinsically fixed - they are either practically 
zero or of 0(1), which is a qualitative difference to mass hierarchies. 

3.2.1. Neutrino mass models 

Starting with mass models, the question to be answered is why the three known ac- 
tive neutrinos have masses of at most 0(1 eV), while even the electron already has 
a mass of 511 keV, and the other SM fermions are by far more massive. In the SM 
itself, neutrinos are even strictly massless, but this is known to be phcnomcnologi- 
cally invalid from the observation of neutrino oscillations, cf. Sec. 3.1. However, the 
reason for neutrinos being massless in the SM is actually a rather artificial one: any 
fermionic mass term would need to couple left-handed (LH) to right-handed (RH) 
fcrmion fields, e.g. C v = —m v VZi>R + h.c, but while (in the 1-family approxima- 
tion) left-handed neutrinos vl appear together with the left-handed electron as 
components of SU(2) doublet fields in the SM, their right-handed counterparts vr 
are simply not included in the particle content. But even introducing right-handed 
neutrinos to create an extended version of the SM does not really cure the prob- 
lem: in that case, we would be forced by the SM-gaugc symmetry to write down a 
SM-like Yukawa coupling term, 

Cy D —LHy v VR + h.c. — s- —L{H)y v VR + h.c. = —TTErnDfR + h.c, (10) 



d Often, the quark sector is disregarded in that context, since quark mixing is not as significant as 
the leptonic mixing. However, even in the more recent literature one can find examples for, e.g., 
the successful prediction of the Cabibbo angle based on symmetry principles. 210-212 
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where L = {v L ,e L ) T is the SM-lepton doublet, H = ia 2 H, H = (H+,H°) is the SM- 
Higgs doublet, and (H) = (0, v) T with v = 174 GeV is its vacuum expectation value 
(VEV) which is generated by electroweak symmetry breaking. Such a mechanism 
induces a mass mrj = y v v for the neutrino, which we have silently given a lower 
index D. This stands for Dirac mass, and it emphasizes that the neutrino in the 
SM-extended by RH neutrinos is a Dirac particle, i.e., distinct from its antiparticle. 
The problem with this way of generating the neutrino mass is that the VEV v is 
fixed to a value of 174 GeV by the known masses of the W- and Z-bosons, and 
unless the Yukawa coupling y v is extremely tiny, y v < 10~ n , the corresponding 
neutrino mass will be too large. 6 Hence, the SM intrinsically predicts a wrong scale 
for the neutrino mass, unless we find an explanation for y u to be tiny. 

The introduction of a Dirac mass reflects the fact that the SM conserves lep- 
ton number at the pcrturbative level. However, lepton number conservation is an 
accidental symmetry of the SM: it has not been imposed, and it is by far not 
sacrosanct. Actually, the SM itself docs violate lepton number, but only at the non- 
perturbative level by so-called sphaleron processes, 213-217 or by higher-dimensional 
operators which are not renormalizable. 218 ' 219 If we accept that lepton number 
may be violated, then gauge symmetry immediately allows a new term in the SM 
extended by RH neutrinos, the so-called Majorana mass term, 

£m - -^{vr) c M r v r + h.c, (11) 

where {vr) c is the charge-conjugate of a right-handed field and, as such, left-handed. 
The factor \ in Eq. (11) arises since the spinors used are Majorana spinors, which 
are effectively real quantities, just as the mass term of a real scalar has a factor 
of i compared to the mass term of a complex scalar. The term in Eq. (11) clearly 
violates lepton number by two units. Even more importantly, the coefficient Mr 
of this term has the dimension of a mass which is not in any way related to the 
electroweak scale v. Hence, Mr can be arbitrary, as long as it is smaller than the 
Planck mass. In particular, Mr will in general not be zero. 

Although there is actually not too much of a reason to suspect Mr to be much 
larger than v, it is often assumed to be. This notion is also the reason for the 
notation Nr instead of vr, which is typically used in the literature and which we 
will also use from now on. If Mr is indeed relatively large, an interesting thing 
happens: the full neutrino mass term can be rewritten in terms of a larger matrix, 

As shown explicitly in Appendix A, this large matrix can be approximately block- 
diagonalized for Mr 3> mrj and leads to a light neutrino mass matrix given by 

m v = —moM^mJ). (13) 



"Note that some authors use the alternative convention (H) = (0,v/\/2) T , in which case v = 
246 GeV. This does, of course, not change the argumentation. 
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What have we gained by this? Since mp <C Mr, a term of O (jf^j is much smaller 
than one, and by this it leads to a suppression of the mass scale mp- If the scale 
of Mr is large enough, say 10 14 GeV, then a mass tud ~ v = 0(100 GeV), as we 
would expect it, would lead to an active neutrino mass of 

™.~^^-«r»Q.v-o.i.v, (14) 

which is just below the upper bounds mentioned in Sec. 3.1. This is the famous 
seesaw type I mechanism. 220 ^ 224 Note that this mechanism can be cast in a Feynman 
diagram, see left panel of Fig. 1: integrating out the heavy neutrino, we exactly 
retrieve Eq. (13). 

This mechanism beautifully exhibits many principles of model building: we try to 
explain certain facts, in this case the smallness of the neutrino mass, by introducing 
new ingredients which are less constrained, in this case the Majorana mass term 
and the heavy fields. At first sight, this might seem to be only a shift of the problem 
to a different sector, or even sweeping it under the carpet. However, the key point 
is that a new model should always lead to certain new predictions, which make it 
testable. In the seesaw type I case, these predictions include the mere existence of 
the right-handed neutrinos, as well as their possible cosmological consequences by 
leptogenesis . 216f By this we can get more insights into the possibilities of Nature: 
many simple-minded models do not work, and our task in model building is to find 
the simplest and most beautiful models which are not only consistent with current 
data, but which also lead to testable predictions, ideally in the near future. It is 
then the task of the experimentalists to measure certain parameters - either more 
precisely or for the first time at all - which we theorists can then use to exclude the 
various candidate models one by one. In the best case, we would by this procedure 
end up with the single model that is as close to the truth as possible. 

To name a few more possibilities, seesaw-type mechanisms can also be achieved 



'Note that, if the seesaw mechanism is associated with an even lower scale, a lot of additional 
cosmological probes might be possible. 225 
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Fig. 1. The Feynman diagrams for the seesaw mechanisms of type I and type II. 
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by adding either a Higgs triplet T (type II 226 ' 227 ) or (typically three) fermion triplets 
£ (type J/J 228 . 229 ^ or D y extending the type I seesaw by additional singlets which 
are not right-handed neutrinos, as in the inverse seesaw. 230 For example, in type II 
seesaw this leads, among other terms, to a triplet Yukawa coupling given by 

Cy T = —Li<T2TyTL c + h.c, (15) 
where the Higgs triplet fields is written in components as 

The point is that the VEV vt is actually induced by the ordinary Higgs doublet 
VEV, due to the term 

Scalar 3 pH T l<J 2 TH + h.C. - M|Tr (T^T) , (17) 

which induces a tadpole term for the Higgs triplet, cf. right panel of Fig. 1. Indeed, 
this leads to a LH neutrino mass as in Eq. (A. 2), where 

_ _ v 2 1_ 

C Yt -> -vlVtVt(vl) c + h.c. = -v L n-j-py T (v L ) c + h.c. = ~-VLm L (v L ) c + h.c. 

(18) 

Note that the potentially large mass Mr suppresses the triplet VEV, which would 
in any case be bound by the deviation of the so-called p-paramctcr, p = M .± CQ ^2 9w 
where M w (M z ) is the mass of the SM W- (Z-) boson, from its SM value of 1 
to vt < 0(1 GeV). 231 If in addition the Yukawa coupling matrix yr has somehow 
small entries, a small LH neutrino mass mi is justified. 

An interesting alternative for neutrino mass models is to suppress the neutrino 
mass not by a tree-level diagram but to let it instead vanish exactly, such that 
it only arises at loop-level. The most simple such extensions of the SM which are 
not excluded are probably Ma's scotogenic model 232 (1-loop) and the Zee-Babu 
model 233 ~ 235 (2-loop). Even higher suppressions are present in the Aoki-Kanemura- 
Seto model 23 ®' 231 (3-loop) or in the Butterfly diagram (4-loop) by Schechter and 
Vallc. 200 ' 201 



3.2.2. The Froggatt- Nielsen mechanism 

What practically all neutrino mass models have in common is that, even though 
they provide an explanation for the smallness of neutrino masses, they usually do 
not give immediate predictions for the mixing angles. 8 Indeed, this is a task that is 
non-trivial to achieve. 

The state of the art of most of the models is to make use of so-called flavor 
symmetries. The idea behind this approach is to extend the undoubtably successful 



s Notable exceptions are, e.g., the Zee-Wolfenstein modei 238 ' 239 or the left-right symmetric exten- 
sion of the scotogenic model. 240 However, the former is excluded by data, and the latter employs 
an additional symmetry, which actually does not make it a pure mass model. 
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application of symmetries, which for example predict the structure of the gauge 
sector of the SM. E.g., the color SU(3)c symmetry in Quantum Chromodynamics 
(QCD) predicts the existence of eight different gluons, and the SU{2)l x U(1)y 
symmetry of the electroweak sector predicts the p-parameter to be exactly one at 
tree-level in the SM. 147 Similar considerations could, in principle, be applied to the 
flavor sector to predict certain mass and mixing patterns. 

The first question to ask is whether the flavor symmetry used should be a con- 
tinuous or a discrete symmetry. Having the SM in mind, a continuous symmetry 
seems to be the natural choice, as such symmetries predict the structures of the SM: 
e.g., the gauge sector exhibits an SU(3)c x SU(2)l x U(1)y symmetry 147 and also 
lepton number could be described by one or more U(l) rotations. 241 Unfortunately, 
a practical problem arises: typically, flavor symmetries are broken for phenomeno- 
logical reasons. One could in principle assume the symmetry to be explicitly broken, 
but then the question might arise in how far the symmetry under consideration is 
actually present in a certain model. Breaking the symmetry spontaneously, how- 
ever, will by the so-called Goldstone theorem 242 ' 243 generically introduce unwanted 
massless scalar particles, which are normally a phcnomcnological disaster. One way 
out would be to gauge the flavor symmetry, as was done, e.g., in Refs. 244, 245. Al- 
ternatively, one could simply use a discrete symmetry, which would not suffer from 
unwanted Goldstoncs at all as long as no accidental continuous symmetry appears 
that gets broken at some point. 

If not the first actual flavor symmetry, at least the most intuitive to understand 
was provided by Froggatt and Nielsen in 1979, by now termed the Froggatt-Nielsen 
(FN) mechnanism. 24e Although originally developed for the quark sector in order 
to simultaneously explain the mass and mixing pattern as well as the occurrence of 
CP violation, 247 the FN mechanism is equally suited for the leptonic sector. The 
basic idea is very simple: since any fermionic mass matrix in the SM arises at the 
electroweak scale v, it is tempting to consider a mass matrix M = vY to originate 
from a Yukawa coupling matrix Y of the form 



where the "natural" entries Y^ at are all of the same order, and the suppression of 
certain elements comes from powers of a suppression factor A, which depend on 
certain generation-dependent "charges" a, and bj. The trick is to re-interpret the 
actual couplings Y^ at X ai+ j as arising from the VEV of a new SM-singlet scalar 
field 6, which is integrated out at a high energy scale A. The "static" couplings are 
hence promoted to dynamic quantities arising from a field which would, in a more 
modern language, be termed a spurion. The charges cij and bj are interpreted as 
being the quantum numbers of a new [7(1 )fn symmetry under which the different 
fcrmion generations are charged differently. Furthermore, the new scalar O, usually 
called flavon h , has a non-trivial charge under ?7(1)fn, and a new set of heavy singlet 




(19) 



h Note that this term is generally used for scalar fields that are charged only under the flavor 
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fermions S a , all suitably charged under the new symmetry, is introduced as well. 

What we gain by the introduction of all these new ingredients is that most of 
the Yukawa couplings, as they appear in the SM, are now forbidden by the J7(1)fn 
symmetry. Instead, new couplings made of higher order terms appear, which couple 
the flavon to either the SM Higgs, one SM fcrmion, and one heavy fermion, or to 
two heavy fermions, e.g.: 

C e , su = Y®I^HSij + h.c. or £ e ,smgiet = + h.c. (20) 

The exact set of terms is determined by the charges chosen. Now, since the flavon 
obtains a VEV (0), this VEV breaks all symmetries to which the flavon couples, 
i.e., the Z7(1)fn- Hence, one can use a "chain" of couplings as the ones in Eq. (20) to 
construct extended seesaw-like diagrams, see Fig. 2 (cf. left panel of Fig. 1), where 
each VEV breaks the Z7(1)fn charge by k units, with k being the charge of the flavon 
field 6. Integrating out each of these heavy fermions introduces powers of the small 
quantity A = and, since these suppressions are different for different SM fermion 
generations, we end up with hierarchical Yukawa couplings at low energies, exactly 
as in Eq. (19). This allows for a first, though simple, understanding of what "flavor" 
actually could be. 

For example, let us take an easy one-generation model of an electron, contained 
in a LH doublet L and a RH singlet cr, where the FN charges of (0, L, e R , H) 
are taken to be (fee, kh, k e , kn) = (1, 2, 0, 0). Then, a term LHe R + h.c. would have 
a C/(1)fn charge of —2 + + = —2 (note the conjugate in the bar!), and hence 
be forbidden. However, if we introduce a new fermion S, we can write down a term 
UlLH^S + h.c. if we take the J7(1)fn charge of S to be k$ = 1 (and also account 
for the correct electrical charge). Equally fine is a term y R SQeR + h.c, which would 
have a FN charge of— 1 + 1 + = 0. The decisive part of the Lagrangian would 
then be 

C D -y L LHjS - VL S^-H*L - y R SOe R - y R eR®*S - MSS, (21) 

where we have written the h. c. -terms explicitly, and we have chosen all couplings to 
be real by absorbing their phases into the fields. We have furthermore introduced a 

symmetry, but are SM singlets otherwise and obtain a VEV. If such singlets do not obtain a VEV, 
they are often termed driving fields or waterfall fields, since they are important for the vacuum 
structure of the scalar potential. 

en S\ Si S3 54 L 
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Fig. 2. Example Feynman diagram for a mass term generated by the FN mechanism. 
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large mass M ol the new fermion S. If now M is so large that S is essentially static 
(i.e., its kinetic term is negligible), then we can integrate out S and S by using the 
Eulcr-Lagrangc equations, e.g., 

I = -»i n - "» e « - MS = d -im <* »• (22 » 



and hence 



Inserting this value back into Eq. (21) and inserting the VEVs (O^*') = w, we 
obtain: 

r 2y L yR T W 2 

C D — LH—e R + h.c. (24) 

If the mass M is taken to be identical to the high scale A, which was not specified 
up to now, then we finally arrive at an electron Yukawa coupling given by 

£ 3 2y L y R X 2 LHAe R + h.c, (25) 

where A = w/A. This has precisely the form anticipated in Eq. (19): the Yukawa 
coupling has a natural size 2yi l y R , which is suppressed by a certain power A i+ e = 
A 2 , due to A being smaller than one. 

Although the basic principle of the FN mechanism is easy to understand, it al- 
ready reveals certain problems that are typically associated with flavor symmetries: 

• Goldstone bosons: 

As mentioned earlier, the breaking the ?7(1)fn m &y introduce undesirable 
masslcss Goldstone bosons. 242 ' 243 To go around that, we either have to 
gauge the continuous symmetry, or otherwise give mass to the Goldstone 
modes by some non-perturbative mechanism as in, e.g., Refs. 248, 249, or 
have the symmetry explicitly broken from the very beginning, as is the case 
for pions. 250 This problem does not only appear when we explicitly impose 
a continuous symmetry from the very beginning, as in the FN case, but also 
if there is an accidental symmetry in the problem. Even worse, we might not 
even be aware of such a symmetry, and wonder why massless states arise at 
all. This problem can be, to some extent, circumvented when using discrete 
flavor symmetries instead, but a model builder should still be careful to 
precisely determine the actual symmetry of the problem. Treatments of 
how to obtain a discrete from the breaking of a continuous symmetry might 
help in that respect, 251-256 and in particular the use of so-called invariant 
polynomials may be of advantage. 257 

• higher-order terms: 

As in Eq. (20), the actual Yukawa couplings at low energies are in fact higher 
order terms. In other words, we are always working in an effective field 
theory (EFT), and hence our Lagrangians arc non-renormalizable. While 
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this may not be a big problem in practice, since one can very well work with 
EFTs (and since the seesaw diagram itself only exists in the framework of an 
EFT), one should still keep in mind that in reality we would need to know 
the fundamental high energy theory in order to truly derive predictions 
such as the sizes of certain couplings from a model. However, since we treat 
the couplings as free parameters, and since not even a plausible candidate 
theory to explain their absolute values is known, we can neglect this problem 
in most practical considerations. 

• unknown high energy sector: 

A bigger problem in practice is the existence of the high energy sector 
itself, in the FN case the flavon (or several flavons) as well as the heavy 
singlet fermions. This is, first of all, an aesthetic problem, since a completely 
unknown sector is introduced but only partially used and often not even 
discussed in what concerns its possible effects or signatures. However, at 
least in the early Universe, this scalar sector could potentially be accessible, 
and by this the associated scalar SM singlet fields could, e.g., be responsible 
for inflation. 258 

• symmetry and charge assignments: 

Already in the FN example, it appears that the charges a, and bj can 
essentially be freely chosen. 1 Not only that, even the whole symmetry group 
can be more or less arbitrarily chosen, and the Z7(1)fn is only an example. 
However, in the spirit of the unification of three generations of fermions, one 
might be tempted to focus on symmetry groups that actually do have three- 
dimensional irreducible representations 3 , as this might signal a common 
origin for, e.g., the three known charged lepton flavors. On the other hand, 
also two-dimensional and even one-dimensional representations are used, 
the latter being applied in the E/(1)fn example. 

The Froggatt-Nielsen mechanism already contains practically all the important fea- 
tures and characteristics of modern flavor symmetry models. In fact, it can serve 
as example throughout the remainder of this review, and actually it is still used 
by many modern models, which are often only variants (though more complicated 
ones) of the "old" idea behind FN. 

3.2.3. Neutrino flavor models 

We will now consider a more modern model, to make clear where we actually need 
group-theoretical aspects. Although the field of flavor model building is an industry 
by itself, we will try to give a flavor of what is going on, which will turn out to be 
very useful in the later sections. A very detailed review on the details of discrete 



'Note that, although typically chosen to have integer values in FN-inspired models, U(l) charges 
can actually be arbitrary real numbers. 
J We shall discuss this term in a second. 
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flavor symmetries by King and Luhn has recently become available in Rcf. 259. 

The model to be presented here is a slight variant the minimal A4 model from 
Ref. 260, which is a particularly economic example. Without aiming to present a 
complete and realistic model, it nicely illustrates how group theory is used to predict 
mass matrices, and hence mixing patterns. 

The basic group theory of A4 is concisely outlined in Ref. 261. The group A4 is 
an alternating group, to be precise, the group of even permutations of four objects. 
Accordingly, it has exactly 4!/2 = 24 2 = 12 elements. All these elements can be 
built out of combinations of powers of two generators, denoted by S and T, which 
fulfill the presentation rules of A4: S 2 = T 3 = (ST) 3 = 1. These rules unambigu- 
ously define the group. It turns out that the group A4 has exactly four irreducible 
representations (irreps), cf. Tab. 1, out of which three are one-dimensional (1, 1', 
1") and one is three-dimensional (3). This essentially means that there are three 
different ways to denote the group elements by (complex) numbers, and one way to 
do that using 3x3 matrices. k Since we want to explain some structure in the mass 
matrices, we make use of the 3 in order to unify the three generations. In other 
words, we simply postulate that the three LH lepton doublets of the SM transform 
under A4 just as the triplet docs, while the right-handed charged leptons transform 
as singlets: 

L = ~ 3, e R ~ 1, m ~ 1", tr ~ 1'. (26) 

Note that this choice is, to some extent, completely ad hoc. An experienced model 
builder would of course know how to assign transformation behaviors in a successful 
way, but we could also take the more naive viewpoint of simply choosing some 
assignments to see where they lead to. 

Table 1. Generators of A4 (with 
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The group theory now comes in because in a Lagrangian we have to combine 
certain fields, in order to obtain the total singlet terms we can include in the La- 



k However, some of these representations are non-faithful: the same number or matrix could be 
used for more than one group element. 



26,2013 3:10 WSPC/INSTRUCTION FILE keV'Merle 



24 Alexander Merle 



grangian. The easiest way to arrive at a neutrino mass is to make use of the (non- 
renormalizablc) dimension-5 Weinberg operator, 21 

A = ^{ViH)za 2 {HL 3 ), (27) 

where A denotes some high energy scale and o~ 2 is the second Pauli matrix. If the 
SM Higgs field H transforms trivially under A4, i.e. H ~ 1, then this operator, 
which leads to a neutrino mass if H obtains a VEV, is in terms of A4 nothing else 
than a direct product of two triplets, 3 ® 3. To know precisely what this product 
is, we need the help of group theory. In the case of A4, it tells us that two triplets 
a = (01, a 2 , a 3 ) and b = [b\, b 2 , b 3 ) are combined as follows, 

3©3 = lffil'ffil"ffi3 s ffi3 A , (28) 

where the resulting singlets are given by, 

1 = (ai&i + a 2 b 3 + a 3 b 2 ), 
1' = (aib 2 + a 2 b x + a 3 b 3 ), 

1" = {a 1 b 3 + a 2 b 2 +a a b 1 ), (29) 

and the two resulting triplets are conveniently decomposed into a symmetric and 
an anti-symmetric expression, 

3s = i(2aifri - a 2 b 3 - a 3 b 2 , -a\b 2 - a 2 b 1 + 2a 3 b 3 , -aib 3 + 2a 2 b 2 - a 3 &i), 

3 A = i(a 2 & 3 - 03^2, aib 2 ~ a 2 b 1: a x b 3 - a 3 b x ). (30) 

Note that this decomposition can always be done, and the two triplets are perfectly 
indistinguishable from an A4 point of view, which we can check easily by acting 
with S3 and T 3 onto a = (a±, a 2 , a 3 ) and b = (bi, b 2 , b 3 ). 

So far, only the trivial singlet combination, the first term in Eq. (29), would be 
allowed in the Lagrangian, which would lead to a relatively boring neutrino mass 
matrix: 

- v 2 — -iv 2 - ( 10 °\ l VeL " 

VOIO/ \v tLj 
(31) 

Even worse, the charged leptons would be completely massless, since any combina- 
tion of an A± triplet a = (cti, a 2 , a 3 ) with a singlet c would generate yet another 
triplet: 

3 (g> 1 = (aic, a 2 c, a 3 c), 
3 (gi l' = (a 3 c, a\c, a 2 c), 

3© 1" = (a 2 c,a 3 c,aic). (32) 

All these combinations would transform non-trivially under A4 and hence cannot 
be allowed terms in a Lagrangian. 
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This problem can be cured by introducing new fields that are again called 
flavons, which are scalar SM singlets but transform non-trivially under the fla- 
vor symmetry. Furthermore, flavons obtain VEVs which can be used to achieve a 
certain structure in the mass matrices. In the model discussed in Rcf. 260, the au- 
thors suggest to use one triplet flavon, 4>s ~ 3, as well as one singlet flavon, u ~ 1, 
which are taken to obtain the following VEVs: 



a A F , (33) 



where Ap is the scale where the discrete flavor symmetry is broken, and the 
parameters as,o are coefficients of 0(1). The triplet flavon field can be used to 
allow for a charged lepton mass term, since 3 <S) 3 = 1 © 1' 1" © cf. Eq. (28), 
and in A4 we furthermore have, 260 

l a ®l b = l (a+b) mod 3 , (34) 

where a, b = 0, 1, 2 and (1°, l 1 , l 2 ) = (1, 1', 1"). For the charged leptons, we obtain 



-> as k e (L e + + L T ) He R + k^ (L e + + L T ) Hfi R + k T (L e + L u + L T ) Ht r 
~- v ' v v ' v v ' 

=1 =1' =1" 

(k e k a k T \ I e R \ 
k e k u k T Lu_r . (35) 
k e k a k T J \T R ) 



Note that the VEV alignment from Eq. (33) makes all contributions look very 
similar; nevertheless, each contribution is a singlet for itself, and hence carries an 
individual coefficient in the Lagrangian. Although the inclusion of the flavon VEV 
leads to a non-zero mass matrix, there will be still two massless fcrmion states, 
since the rank of the mass matrix in Eq. (35) is one. This suggests that actually one 
would like to introduce yet another flavon to give masses to the charged leptons, 
just as done in Rcf. 260 - see that reference for more details. 

For the neutrinos, in turn, we can make use of both flavons. Since a Majorana 
mass matrix as the one originating from Eq. (27) must be symmetric, cf. Appendix 
A.l, the only triplet to combine with the flavon 4>s is the symmetric one. Applying 
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Eq. (30) yields 

1 — ,TTT\(<fe) _ a sv 
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Similarly, the trivial singlet combination in Eq. (27) can be married with the singlet 
flavon u, yielding 
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= ^K c i; ^,^)|00l| |^| . (37) 

Making the redefinition a.g/3 -> ag, we finally arrive at 

1 _ (u) v 2 /«o + 2a s 

x (Lcfl)(HL)g D X (^ C L ,^,^ L ) -as 

(38) 

We have now seen how the structure of the mass matrices in Eqs. (35) and (38) 
arose. The next step would be to diagonalizc both mass matrices, and to determine 
the PMNS-matrix by the mismatch of the two diagonalizations. A very detailed 
treatment of this point can be found in Refs. 262, 263, and some more information 
is also given in Appendix A. 2. 

Let us finish by noting that the neutrino mass matrix in Eq. (38) is form- 
diagonalizable, i.e., it can always be diagonalized by the same matrix U u , no matter 
which values the parameters as,o have. As it turns out in the particular case at 
hand, the matrix that does this job is tri-bimaximal, 264 which would lead to mixing 
angles Q\ 2 = arctan {^^j i $13 = 0, #23 = f • Since a vanishing angle $13 is, however, 
excluded, cf. Eqs. (4) and (5), we will need a non-trivial charged lcpton mixing, 
U e 7^ 1, in order to obtain a phenomenologically viable mixing matrix C/pmns = 
XJ\\J V . This would, again, make it necessary to modify the charged lepton mass 
matrix, by introducing even more fields that ideally should not destroy the form of 
the neutrino mass matrix. The reader should at latest now appreciate that finding 
a working model is indeed a non-trivial task. 

Note finally that in our calculation we have relied heavily on the chosen vacuum 
alignment, cf. Eq. (33). However, we have not written down a full scalar poten- 
tial including all allowed terms involving the SM-Higgs as well as both flavon fields. 
Even if we succeed in writing down such a potential, it is a highly non-trivial task to 
show that Eq. (33) is at all a minimum of this potential (showing the even stronger 
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condition that it is actually the global minimum is a task that is practically im- 
possible to tackle in many realistic models). The typical way to proceed is to do 
a numerical study in order to show that the alignment chosen is at least fine for 
a certain choice of parameters. We will here not enter such details, but we should 
keep in mind that in a realistic model we would actually need to do that. 

We will now see how the concepts we have discussed come into play when trying 
to find working models that include kcV sterile neutrinos. Not only will we have 
to convincingly explain all neutrino data, but we will also have to be careful to 
be in agreement with certain bounds that are specific for keV neutrinos, such as 
the non-observation of the astrophysical X-ray line (cf. Sees. 2.3 and 4.1), and in 
particular the generation of the correct DM abundance. 

4. General features of settings with keV Sterile Neutrinos 

In this chapter, we will summarize a few theorems, conditions, and other features of 
working models for keV sterile neutrinos. This collection will help us to appreciate 
many of the steps which are implicitly undertaken in concrete models. To mention 
one example, a typical argument used against keV sterile neutrino models is that 
one should be skeptical if the seesaw mechanism (cf. Appendix A) actually works, 
since after all one would have to divide by a small mass of O(keV) when applying 
the seesaw formula, Eq. (13). While this is certainly a dangerous enterprise at first 
sight, one can actually show that in any working model of kcV sterile neutrinos which 
respects the X-ray bound, this is no problem at all. A couple of such properties, to 
be respected by a large class of models, will be presented in the following. 

4.1. The X-ray bound 

Although already mentioned in Sec. 2.3, we will now stress the importance of the 
X-ray bound, arising from the decay N — > ^ a is a j and the non-observation of 
the corresponding astrophysical line. Note that we have now generalized the keV 
sterile neutrino N to originate from any generation i, and that we have furthermore 
indicated that any active neutrino flavor a = e,u,r could be produced. However, 
since we can never observe this flavor we have to sum over all possibilities. 

In order to get a better understanding of the constraint arising from this decay, 
we have depicted the two lowest order 1-loop Feynman diagrams for the reaction 
Ni — > v a l in Fig. 3. As already mentioned in Sec. 2.3, the diagrams look very 
similar to a fi — > e7 transition, which is even discussed in some textbooks (see, 
e.g., Ref. 265). More formally, the process is a special case of the general radiative 
fcrmionic decay fi — > /27- 104 Since the charged lepton and the VF-boson inside the 
loop both are electrically charged, the external photon can couple to either of them, 
which is why we have actually two diagrams for each active neutrino flavor v a . The 
decisive point in the diagrams is the left vertex: the mass eigenstate sterile neutrino 
Ni is not purely right-handed, but it has a small left-handed active admixture, 
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which is parametrized by the (small) active-sterile mixing angle 9 a i, cf. Appendix 
A. 2. Hence, there is a possibility for this mass eigenstate to couple to a charged 
lepton e a and a VF-boson, but this coupling is suppressed by a factor 9 a i and the 
amplitude M. a i for Ni — > v a ^j is proportional to exactly the same factor, as indicated 
in the diagrams. 

Note that we can produce any active neutrino flavor a with a probability pro- 
portional to |.M Q i| 2 . However, since we have no way to measure the final state 
flavor a, we have to sum over all possibilities, a — e,/i, r. Note that, although we 
cannot determine the flavor of the final state active neutrino, we nevertheless know 
by laboratory experiments that the different flavors v e ,^,T are fundamentally dis- 
tinct. This is why we have to sum over probabilities {incoherent summation) rather 
than over amplitudes [coherent summation), differently from what we did for, e.g., 
neutrino-less double beta decay, cf. Sec. 3.1.2. Thus, the decay rate Ni — > into 
all possible active neutrino flavors must be proportional to the following quantity, 



which is called the i-th active sterile mixing angle square and usually denoted as 
&f . This is the quantity we can put an upper bound on by a non-observation of the 
monoenergetic photon 7. 

The precise bound originating from different satellite experiments can be found 
in Refs. 29, 73 (based on Refs. 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 
115, 116). For our purpose, a simplified version of this bound as used in Ref. 24 is 
perfectly sufficient: 



Note, however, that more recent non-observations of the X-ray line for certain galax- 
ies yield even stronger bounds, cf. Refs. 117, 118. Corresponding updates of the 
simplified bound in Eq. (40) are available. 266 

Let us now get a more precise understanding of the connection between the 
active-sterile mixing and the entries in the full neutrino mass matrix. As we have 
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Fig. 3. The Feynman diagrams for the radiative decay of the sterile neutrino, Ni — > v a ~f. 
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just seen, the definition of the i-th active-sterile mixing angle 9i is 



\e 



where 9 m = U a ^ +l = 



i* D M R > 



V R 



(41) 



Note that we have expressed the generation-dependent active-sterile mixing 9 a i in 
terms of the full neutrino mixing matrix U as defined in Appendix A. 2. In the basis 
where the RH neutrino mass matrix is diagonal (and real), Mr = diag(Mi, M%, M 3 ), 
we have Vr = 1, and the above formula simplifies to 
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This can be simplified further, 
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We can now see how to understand the active-sterile mixing 8 2 : it is simply the 
proportionality factor between the sterile neutrino mass Mj of JVj and the corre- 
sponding active-neutrino mass rrii of Vi in a seesaw type I framework, cf. Appendix 
A. 2. Hence, it should also be intuitively clear that 8 2 indeed parametrizes the small 
"active fraction" of the otherwise sterile mass eigenstate Ni , and it can thus be used 
to quantify how often a neutrino which is mainly sterile can decay via a coupling 
to active partners, such as the charged leptons and the W^-bosons in Fig. 3. 



4.2. Seesaw theorems 

Proceeding with some more technicalities, we will now shortly discuss two very 
useful theorems which hold for seesaw type I situations. The first theorem 267 deals 
with the potential issue of applying the seesaw formula, Eq. (13), to models involving 
keV sterile neutrinos: 

Theorem 4.1. keV seesaw practicality theorem 267 

For any working keV neutrino model (i.e., it is consistent with the X-ray bound), 
the seesaw formula works. This conclusion can only be altered if the decay mode 
N — > 1/7 does not exist (or is not effective) for some reason. 

Proof: The active-sterile mixing angle 9 a i is given in Eq. (41). In a basis where Mr 
is diagonal we have 8 m = (m j^"' , cf. Eq. (42). If M 1 = O(kcV) < M 2 , 3 , then the 
only critical property to be fulfilled for the seesaw to work is \8 a i\ <C 1. From the 
rough X-ray bound, Eq. (40), in combination with the rough Lyman-a bound, Mi > 

1 keV (cf. Sec. 2.1), it follows immediately: (m^mD)^ < 1.8- 10~ 5 keV 2 (tj^) < 

1.8-10" 5 keV 2 < Mi, which implies \m Del \ 2 + \m D + \m Dr i\ 2 < M 2 , and hence, 
due to the sum over absolute values, \mD a i \ ^ Mi for a = e, \i,r. □ 
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Thus, as long as a model for keV sterile neutrinos respects the X-ray bound, one 
can always apply the seesaw formula. If a certain model violates the X-ray bound 
it is either excluded, in which case we do not need to talk about it, or the decay 
mode Ni — > v"f is somehow forbidden or strongly suppressed, e.g. by a symmetry 
stabilizing the keV neutrino N±, cf. Ref. 268. In the latter case one might have to 
check the validity of the seesaw formula, if applicable, but in most models one can 
blindly use it if the X-ray bound is known to make no problems. Note that the 
above argumentation is not changed if there are additional type II contributions 
to the light neutrino mass matrix, since any terms in mi do not contribute to the 
active-sterile mixing, cf. Appendix A. 2. 

The second theorem was first mentioned in Ref. 269 and later on brought into 
a more modern form in Ref. 270. It is very useful for checking the consistency of 
diagonalizations of complicated neutrino mass matrices: 

Theorem 4.2. Seesaw /airplay rule 269. 

In a seesaw type I setting with p LH doublets and q < p RH neutrinos, one obtains 
at most q light neutrinos with a non-zero mass, while at least {p — q) neutrinos 
remain massless. In particular, in a setting with 3 LH doublets and only 2 RH 
neutrinos, one light neutrinos is exactly massless. 

Proof: The proof of this theorem is trivial if one visualizes the structure of the 
corresponding matrix. In Eq. (12), the upper left zero block is apxp matrix. In order 
for the first p column vectors of the full mass matrix to be linearly independent, 
one would need q > p components to form enough unit vectors. Due to q < p, 
this is never possible and at most q of the first p column vectors can be linearly 
independent. Even if the remaining q column vectors consisting of the right mo 
and Mr blocks are all linearly independent, the maximum rank of the full mass 
matrix is 2q. Hence, the matrix has at most 2q non-zero mass eigenvalues, and at 
least (p — q) of the light neutrinos remain massless. □ 

This theorem is exactly what is at work in, e.g., the model to be presented in 
Sec. 5.2.2: since the keV sterile neutrino is decoupled from the mass matrix, the 
lightest neutrino remains exactly massless in that approximation, cf. Eq. (59). Note 
that, however, a small but non-zero mass is nevertheless generated at 2-loop level 
by diagrams involving only SM-fields. 271 

4.3. Cancellations in neutrino-less double beta decay 

We also want to present a theorem discovered in Refs. 272, 273, which reveals an 
interesting effect of light sterile neutrinos on the effective mass \m ee \ in neutrino- less 
double beta decay, cf. Eq. (7). 



Theorem 4.3. Ov/3/3 cancellation theorem ' 73 

In a seesaw type I situation where all RH neutrinos have mass eigenvalues below 
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100 McV, the effective neutrino mass \m ee \ measured in neutrino- less double beta 
decay vanishes exactly. If the neutrino mass matrix is in addition form dominant, 
then this cancellation happens piecewise for each generation. In such a case, even 
having only a few RH neutrino masses below 100 MeV leads to partial cancellations 



Proof: Whenever the mass of a neutrino is below the nuclear momentum 
transfer in OV/3/3, \q\ ~ 100 McV, it contributes to the decay amplitude by light 
neutrino exchange, 177 and this contribution is given by U ek mk- Hence, for three 
LH and three RH neutrinos with all masses below \q\, the effective mass is given by 

\ m ee\ = ELi u lk m k + ELi U e,3+k M k = [ M ^ x6 ] ee , whicn is zero by definition 
in a type I seesaw setting. If the mass matrix is in addition form dominant, cf. Ap- 
pendix A. 4, then Eq. (A. 35) holds in the basis where Mr is diagonal. One can imme- 
diately conclude that (m* D ) e k = i\J m k M k XJ ek and that U e .3+k = &ek = 

cf. Eq. (41). One readily obtains ul 3+k M k = 8 2 ek M k = -^U* k M k = -U 2 ek m k . □ 

Note that this theorem implies in particular that, in models with one keV sterile 
neutrino N\ and two other neutrinos A*2,3 with M 2 ,3 > \q\ — 100 MeV, the effective 
mass for a form dominant setting is given by 266 

|m ee | = \m 2 s 2 12 cl 3 + m 3 s 2 3 e^- Q) | , (44) 

and similar if the keV sterile neutrino originates from a different generation. This 
cancellation was not recognized in Refs. 274, 275, which contained the first discus- 
sion of 0vftf3 in the context of the i^MSM. 1 Instead, it was shown that the contribu- 
tion from the keV sterile neutrino is tiny, which is why it was neglected. Hence, the 
cancellation could not have possibly been taken into account in these references. In 
particular, the formula for \m ee \ obtained in in Rcf. 274 for the case of inverted mass 
ordering would be wrong if indeed Ni was the keV neutrino and form dominance 
held: instead of the standard contribution, as obtained in Ref. 191, the correct ex- 
pression would be given by |m ee | ~ -\/| Am^ |s 2 2 c 2 3 , which is obtained from Eq. (44) 
for the limit m 3 < mi < 1112 with 7713 — > 0. 



4.4. A little No- Go theorem 

The final theorem has appeared in a slightly simplified version in Refs. 23, 27 and 
was later on generalized in Ref. 267. Note that in the theorem and in the proof we 
will make use of the so-called Casas-Ibarra parametrization, the details of which are 
explained in Appendix A. 3. 

Theorem 4.4. No- Go theorem 23 ' 27 ' 267 

A quasi-degenerate light neutrino spectrum and keV sterile neutrino Dark Matter 



'in particular, this references focused on the "standard" production mechanisms (DW and SF) of 
keV sterile neutrino DM, in which case the lower bounds on the keV mass are strongest. 
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contradict each other in a seesaw type I setting with a real Casas-Ibarra matrix R, 
if the sterile neutrino mainly decays via N% — > v^j. 

Proof: Inverting the known Casas-Ibarra formula, Eq. (A. 22), one obtains 
mi) = iJ7*diag( v / mT, y/rn^, v /m^)R T diag(y/Mi, y/M 2 , y/M 3 ). Using Eq. (42) 
yields 6 m = ~iY.k,L m u c<k^/m^8 kl R\ m ^=8 mi = ypjjfcUckdkiRu = 

- i Y.k\[W z u *kR* lk - With the hel P of E q- ( 41 )> one can show that d i = 

Eq l^ml 2 = Eq E/U ^M~ UakU al R ik R H = Efc,Z ^ A//"' (^q U L U akj R il R ki = 

Efc Tf R ik R \i- For quasi-degenerate light neutrinos, one has rrik — to for k = 
1,2,3, and if the orthogonal matrix R is real then R^ = R T and Efc -^-ifc-^L = 
&u = 1. For i = 1, one finally arrives at Q\ = -^2-. Applying the X-ray bound, 
Eq. (40) yields a relation that cannot be fulfilled: applying the lower limits on M\ 
from structure formation requires too <^ Am\ [cf. Sec. 2.1 and Eq. (40)], while 
quasi-dcgeneracy would require too 3> \J Am A . □ 

The detailed implications of and ways around the No-Go theorem are discussed 
in Rcf. 267. The important point we want to stress is that R is real in particular if 
CP invariance holds. Hence, models for keV sterile neutrinos with CP invariancc 
will be in trouble if they predict a quasi-degenerate light neutrino spectrum. On 
the other hand, if the light neutrinos were experimentally known to have a quasi- 
degenerate mass pattern and if at the same time keV sterile neutrinos were the Dark 
Matter, this would implicitly prove the existence of CP violation in the neutrino 
sector, as long as the light neutrino mass is generated by a seesaw type I mechanism. 

Equipped with these general considerations and theorems, we are now ready 
to enter the central discussion of this review, exploring the known models for keV 
sterile neutrino Dark Matter. 

5. Models for keV Sterile Neutrinos 

We will now turn to the core chapter of this review, where we will give a supposed- 
to-be-complete discussion of the models of keV neutrinos on the market. Hereby, 
we will not follow the historical order in which the different models have been 
developed, but we will rather apply a more pedagogical way through the jungle of 
possibilities, in order to give the less experienced reader a solid guideline, and to 
provide the more experienced reader with a maybe slightly different viewpoint. 

In the course of the chapter, we will strongly distinguish in terminology between 
models and scenarios. In the definition used here, we will call a certain setting 
scenario, whenever it can accommodate for a keV sterile neutrino, but does not 
give any explanation for the appearance of the keV scale. One typical example for a 
scenario is the j/MSM, 23 although it is actually termed "model" in the corresponding 
reference. On the other hand, we call a certain setting a model whenever there is 
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an explanation for the appearance of the keV scale or, rather, for a suitable mass 
hierarchy or the existence of a suitable new scale. One example of what we would 
call a model would be a L e — — L T flavor symmetry extension of the ^MSM, as 
presented in Rcfs. 276, 277, 278. 

The reason for this distinction in terminology is motivated by the usage of the 
terms in contemporary elementary particle physics, however, admittedly the terms 
can overlap in certain cases. We nevertheless try to stick to this terminology as 
strictly as possible, in order to give better guidance to less experienced readers. 
But we also stress that the distinction in terminology does not involve any kind 
of judgement. Indeed, both concepts are extremely useful in their respective con- 
texts: scenarios, on the one hand, are very useful for phenomcnological studies, as 
the recent very complete treatment of the phenomenology of the i/MSM clearly 
illustrates. 29 Models, on the other hand, add the necessary building block of an ex- 
planation for the existence of the keV scale (or, rather, of a suitable scale hierarchy), 
which is important since the keV scale is to our current knowledge unrelated to any 
other fundamental scale in physics. Such models typically have "back reactions" on 
several phenomcnological parameters, as e.g. explained in Rcf. 279. 

Naturally, the ultimate goal we should work towards is to combine scenarios 
and models, in order to identify the most promising theory for keV sterile neutrino 
Dark Matter. This "ultimate" theory should involve a convincing production mech- 
anism as well as a proper explanation for the appearance of the keV scale, and its 
phenomenology should agree with all cosmological observations and particle physics 
experiments. On top of that, any candidate theory should yield testable predictions 
which can be probed in terrestrial and/or extra-terrestrial experiments. While e.g. 
the r/MSM does not have too many direct tests, certain models manage to entangle 
the existence of the keV scale with solid predictions for low energy neutrino ob- 
servables such as mixing angles or effective masses. These models offer testability, 
not only by their phenomenological predictions but already when trying to combine 
them with certain production mechanisms. 

5.1. General mass schemes for keV neutrinos 

Before discussing actual models, we would first like to illustrate the two main ap- 
proaches to explain the appearance of the keV scale. As we have already mentioned, 
a mass of a few keV is, to the best of our current knowledge, not connected to any 
"fundamental" energy scale in Nature, such as the Planck scale or the electrowcak 
scale. This is a pity, because in a quantum field theory it is intrinsically difficult 
to predict absolute scales. What we can predict, however, are hierarchies of scales, 
i.e., why a certain scale is larger or smaller than another one. While some of these 
explanations may sound a bit artificial, this is nevertheless the best way we know 
to predict scales at all. Furthermore, it is a well-known principle to relate "new" 
energy scales to known ones, just as the energy levels of an atom are intrinsically 
connected to the size of the Rydbcrg energy. 
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The two main mass shifting schemes that are, in various versions, generically 
used to "explain" (or at least "motivate" ) the existence of a keV scale are depicted 
in Fig. 4. For obvious reasons, we will call these two main schemes the Bottom-up 
scheme and the Top-down scheme. These two schemes nicely illustrate the two main 
ways to arrive at the "in-between" scale of a few keV: either the mass of the sterile 
neutrino that is supposed to be the Dark Matter is actually zero, but this natural 
value is corrected by some model-specific mechanism to yield a non-zero mass of 
O(kcV). Or the fundamental mass Mr of the RH-ncutrinos is actually much higher, 
and for some reason the mass of one sterile neutrino is suppressed to yield a physical 
mass of only a couple of keV. 

Note that we have, in passing, used the terms right-handed neutrino and sterile 
neutrino to be practically equivalent. This is okay, as long as we are talking about 
SM charges only: the right-handed neutrino is a total singlet under the SM gauge 
group, and hence it is sterile. However, we have to keep in mind that there are two 
loopholes in this terminology: 

• First, the physical particle (which is a mass eigenstate!) will actually be 
neither purely right-handed nor purely sterile, due to the structure of the 
mass terms and the active-sterile mixing, cf. Appendix A. 2. However, the 
standard terminology is to refer to the mass eigenstates as sterile neutrinos 
and to the chiral fields in the Lagrangian as right-handed neutrinos. Nev- 
ertheless, both terms are often used in a more or less equivalent manner 
in the literature, and in many cases one has to conclude from the context 
which physical meaning is actually referred to. 

• Second, even though the term sterile is used, this only refers to SM- 
interactions. As soon as we go beyond the SM by extending the gauge 
group, the RH-neutrinos will not be total singlets anymore, in general. For 
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Fig. 4. The two generic mass shifting schemes for keV sterile neutrinos, in a setting with three 
right-handed neutrinos. Typically, Ni is taken to be the keV sterile neutrino, but there are models 
where this is not true. 
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example, in left-right symmetric models, 88,280 where the SU(2)l x U(1)y 
symmetry of the SM is extended to SU{2) L x SU(2) R x U(1) B -l, right- 
handed neutrinos arc charged non-trivially under the SU{2)r x U(1)b-l 
subgroup. Furthermore, even if the gauge group is not extended, the RH 
neutrinos can easily have relatively strong interactions with non-SM scalar 
particles, as for example in the scotogenic model. 232 

Nevertheless, as done in the bulk of the literature, we will use both terms in a 
pretty analogous manner, not to hinder the flow of the text. Keeping the above two 
remarks in mind is perfectly sufficient for the purpose of understanding the models. 

In addition we want to remark that in the schemes in Fig. 4 we have assumed the 
existence of three right-handed neutrinos ^V 12 ,3, which have masses M li2 ,3- Out of 
these, our goal is to find an explanation for why M\ (or whichever is the keV sterile 
neutrino mass) has a value of a few keV. Then, the corresponding mass eigenstate 
field N± (or A^) is referred to as the ke V sterile neutrino. However, in the literature 
one can also find examples with more or less than three right-handed neutrinos (see 
Ref. 54 for an exhaustive collection). While in such cases, the schemes in Fig. 4 
would of course have to be altered (by adding or removing some of the energy 
levels), the basic principles illustrated in the two schemes nevertheless remain the 
same: in practically all cases known to the author, either a zero mass is lifted or 
a larger mass is suppressed, unless no completely new scale unrelated to anything 
else is introduced and then more or less arbitrarily set to O(keV). 

Equipped with the imagination of the schemes, as well as the terminology to be 
used, we are now prepared to enter the zoo of the actual models, which all attempt 
to find an explanation for the existence of the keV scale. 

5.2. Models involving the Frog gatt- Niels en mechanism 

First we will discuss models involving the Froggatt-Nielsen mechanism. There exist 
pure FN models, which use only FN, and mixed FN models, which employ the FN 
mechanism in combination with, e.g., a discrete flavor symmetry. 

5.2.1. Pure FN models 

A very complete discussion on how to build models for keV sterile neutrinos using 
the FN mechanism was given in Ref. 281. As we have already seen in Sec. 3.2.2, 
the FN mechanism can modify a mass matrix or Yukawa coupling matrix by pow- 
ers of a small number A. Hence, we can immediately conclude that FN inspired 
models resemble a typical top-down setting: the "natural" values of the masses re- 
ceive generation-dependent suppressions, which lead to small values of the physical 
masses, as illustrated in Fig. 5. 

The relatively general approach discussed in Ref. 281 uses the same ingredients 
as the models presented in Rcfs. 282, 283. The minimal particle content to arrive 
at a model with the power of analytical predictions is given by the SM leptons, 
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together with three right-handed neutrinos iVi.2,3 and two flavons Oi,2, which are 
charged under the J7(1)fn symmetry as well as an auxiliary Z 2 symmetry: 111 

61,2: {01,03-,+,-), 
-^1,2,3 : (/lj /2, /3i + , +! — )j 

eTa3 ■ (ki,k 2 ,ks; +,+,-), 



#1,2,3 : (51,52,53;+,+,-), (45) 

The most general (i.e., seesaw type II) Lagrangian that leads to masses in the lepton 
sector is then given by 

a+b=ki+fj a+b=gi+fj 

£ = - yi j emHL jL X^X h 2 + h.c. - J2 Y% N^H L jL XlX b 2 + h.c. (46) 

a.b,i.j a.b,i.j 
a+b=fi + fj a+b=gi+gj 

- E ^L) E m i iL jL X a 1 X b 2 + h.c. - ^(N^y c M l iN jR X^ + h.c., 

where all Yukawa and mass matrix entries denote their natural values, while all 
suppressions come from FN factors. These suppression factors depend on the flavon 
VEVs (0j) as well as on a high energy scale A, cf. Eq. (25), 

A. - <f>. (47) 

Although the flavon VEVs can be complex, all but one phase can be rotated away. 282 
Hence, it is sufficient to use one real parameter A and one complex parameter R, 

A =T' "-fif-**-- <48) 



11 As explained in Ref. 282, this auxiliary symmetry is necessary to have a model with CP violation. 
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Fig. 5. Mass shifting scheme for the FN mechanism. The generation-dependent FN-charges g; 
strongly suppress the lightest sterile neutrino mass Mi, which clearly resembles the top-down 
scheme, cf. right panel of Fig. 4. (Figure similar to Fig. 2 in Ref. 281 ) 
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where Rq and cxq are real numbers. Note that ao is responsible for all (Dirac and 
Majorana) CP violation in the model. Using this parametrization, we get 

^ E (x)°(T) i = AO+W ' (49) 

The goal of any model is to reproduce the known data, in our case the charged 
lepton masses. A successful choice of FN charges (with a = 0, 1) is given by 282,284 

01,2 : (-1,-1), 
£1,2,3 = («+ l,a,a), 

eT^i: (3,2,0). (50) 

The important point is that these assignments do not constrain the FN charges 
of the RH-neutrinos, (31,32,33)- Indeed, it is trivial (though a bit cumbersome 
without the use of a computer algebra package) to show that the light neutrino mass 
matrix arising from the Lagrangian in Eq. (47) is independent of gi. The explicit 
proof for the case of only one flavon field is given in Ref. 281. As explained in 
Refs. 282, 285, the reason is simply that the seesaw mass term, due to its Majorana 
nature, breaks any global symmetry of the RH-ncutrino fields. In particular this is 
true for any assignment of a lepton number, which is typically some kind of global 
U(l) symmetry 241 and causes the corresponding mass term to be lepton number 
violating. In the same way, however, the J7(1)fn is broken by that term, and hence 
the RH-ncutrino charges (31,32,53) cannot appear in the seesaw formula. 

It is exactly this fact that was exploited in Ref. 281: the charges (51, 32, 33) can 
be freely chosen to generate a suitable RH-neutrino mass pattern, and the seesaw 
formula is always guaranteed to work, in spite of the presence of potentially very 
light RH-neutrinos whose masses are "divided by" in the seesaw formula. This is 
the key point why FN-inspired models are extremely well suited to explain low-scale 
sterile neutrinos, and in particular keV sterile neutrinos. 

In the j/MSM, 23 which is more or less the minimal scenario for keV sterile neutri- 
nos, the "heavy" sterile neutrinos N23 need to have masses of a few GeV. Depending 
on the scenario under consideration, one might also be interested in higher masses. 
In any case, for the scenarios under consideration, a suitable mass pattern may arise 
for the second and third generation charges being larger than g\ by at least three 
units, gi > 32,3 + 3. n Having this requirement in mind, two "minimal" scenarios were 
identified in Ref. 281: Scenario A where (31,32,33) = (3,0,0) and Scenario B where 
(31, 32, 33) = (4, 1, 0). Denoting the RH mass scale by M and assuming the natural 
mass matrices to be democratic, = Mq Vi,j, these scenarios immediately lead 



"Note that this statement may be changed if one attempts to find models including sterile neutrinos 
at the eV scale. Such examples are discussed in, e.g., Refs. 279, 286. This is one instance where 
the specific requirements for models with eV or keV sterile neutrinos can exhibit considerable 
differences. 
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to the following sterile neutrino mass patterns: 



A(3, 0, 0) : Mi = M A 6 2i2§y 1 + R% + 2i?g cos(2a ) , 
M 2 = M , 

M 3 = M (1 + A 6 [l + i? 2 (3cos(2a ) + 3i?g cos(4a ) + R% cos(6a )]) , 
B(4, 1,0) : Mi = M A 8 2R^ I + i? 8 - 2R% cos(4a ) , 
M 2 = M A 2 , 

M 3 = M (1 + i?o A2 cos(2a )) . (51) 

Indeed, these two scenarios reveal the desired structure: if A ~ 0.1, then the require- 
ment Mi = 0(1 keV) leads to M ~ (10 6 kcV - 1 GcV, 10 8 kcV - 100 GcV) for 
scenario (A,B). Even larger values of the FN charges would easily be able to push 
Mo to even higher values, due to the exponential dependence on the FN charges: 
already for g\ = 5, one would achieve Mo ~ 10 TeV. 

The two scenarios (A,B) are combined with the two possible charge assignments 
for the LH lepton doublets, Assignment 1 (a = 0) and Assignment 2 (a = 1), cf. 
Eq. (50), which lead to two different possibilities for the charged lepton mass matrix, 

/y ( n i? 2 ,4A 3 - 4 y e 12 i? 2 A 2 - 3 y e 13 5 , 2 i*A 2 - 3 \ 

Mi 1 ' 2 ) = v y e 21 B 2 A 2 < 3 Y 22 B 02 \ 1 > 2 F e 23 i?A 1 ' 2 . (52) 

Xy^rx 1 ' 2 o, r e 32 i?A f^a ' 1 / 

Again assuming the natural mass matrices to be democratic, Y£ 3 = Y e Vi,j, one 
obtains e.g. for Assignment 1 the charged lepton masses 

m e = mo A 3 Rq , 



■in 



ii 



= m A + A 2 Rl cos(2a ) + fl "~f" +3 \ , (53) 



m T = m (1 + f-R 2 ^ 2 ) I 

where mo = vY e . Using the measured values for m e , m M , and m T , one can calculate 
the mass ratios to obtain A ~ 0.06 and Rq ~ 1.18 at lowest order in A. The phase «o 
is not constrained to that order, but choosing ao = 0.67 causes the 0(A 3 ) correction 
to m M /m T to vanish exactly. As typical for seesaw models, 28 the optimal choice for 
A turns out to be a bit smaller than the standard choice of 0.22. 288 

As visible in Eqs. (51) or (53), the FN mechanism is extremely well suited to 
predict the desired mass hierarchy. However, the FN symmetry is a U(l) group which 
has only l-dimcnsional representations. Unfortunately, these are not sufficient to 
predict more eleborate structures of the mass matrix as, e.g., forcing certain entries 
to be exactly equal. One way out is to depart from the assumption of democracy, 
as illustrated in Ref. 281. In practice, the most convenient way is to simply write 
numerical coefficients of 0(1) in front of all Yukawa and mass matrix elements. 
These coefficients are not predicted by the FN mechanism, which can be seen as 
disadvantage. On the other hand, the FN mechanism nevertheless imposes a certain 
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structure on the mass matrices, which makes it actually very easy to find suitable 
coefficients that lead to consistency with all experimental values. 

To give one explicit example, we want to shortly illustrate Model 1AI from 
Ref. 281, which is based on Assignment 1, Scenario A, and a seesaw type I mech- 
anism. By imposing the measured values of the charged lepton masses, of the neu- 
trino mass square differences, as well as the previously determined values of the 
parameters (A, Ro,ao), the numerical coefficients have been determined such that 
mixing angles within their (at that time) 3er ranges were predicted. In the case of 
Model 1AI, the mass matrices found are given by: 

/Y c n B 2 X 3 Y C } 2 B 2 \ 2 Y c } 3 B a RX 2 \ /0.81B 2 A 3 l.UB 2 X 2 0.29RX 2 ' 

AfW = v Y 21 B 2 X 2 Y 22 B Q X Y 23 RX -> M e0 2.005 2 A 2 1.13A 2.50i?A 
\Y 31 RX F e 33 / \3.71i?A 0.35 

/M^B e X 6 MfB 2 X 3 M 13 RB 2 X 3 \ /0.38B 6 A 6 O.3IB2A 3 \.2&B 2 RX 3 

il4 A) = • MfB -> M 0.31B 2 A 3 4.18 

V • • Mf ) \l.26B 2 RX 3 4.81 

/Y^BA 4 Y^ 2 B 2 X 3 Y^ 3 RB 2 X 3 \ /0.75B 4 X 4 0.15B 2 X 3 1A2B 2 RX 3 

m% A) =v\y^B X Yg 2 -> m D o 0.51A 0.13 

XYj^RX Y^ 3 J \3.32i?A 2.93 

where B 2n = 1 + R 2 + ... + R 2n , and M e o, Mq, and moo denote the characteristic 
mass scales of the charged leptons, RH masses, and Dirac masses, respectively. Note 
that some entries are zero, since they are forbidden by the auxiliary ^-symmetry. 

While Model 1AI still does not predict the absolute light neutrino mass scale 
m v , it is powerful enough to predict their mass ratios, their mass ordering, as 
well as all leptonic mixing angles and phases. Calculating all these quantities from 
the mass matrices of the model, one obtains: 281 sin 2 9\ 2 = 0.28, sin 2 613 = 0.018, 
sin 2 6> 23 = 0.54 , 5 = 4.21, a = 0.58, (3 = 1.25, m 1 /m v = 0.0014, m 2 /m v = 0.19, and 
ms/rriv = 1.03, which means normal mass ordering. 

Note that, although the FN mechanism is not very well suited to predict more 
complicated leptonic mixing patterns, it nevertheless intrinsically predicts that 
$13 7^ in the absence of extremely unlikely cancellations. This is in contrast to 
some discrete flavor symmetries with a tendency to predict vanishing #13. This is 
particularly remarkable since the predictions from Ref. 281 appeared before the ac- 
tual experimental determination of #13. Nevertheless, even today, Model 1AI is valid 
within the 3er ranges of the leptonic mixing parameters. 

Let us end this section by mentioning an interesting feature of the FN mecha- 
nism, which has been discussed in Ref. 281. While the C/(1)fn charge assignments 
seem quite arbitrary at first sight, there nevertheless exist many settings that do not 
work in connection with FN, in particular in the context of keV sterile neutrinos. 
Out of all the requirements mentioned in Ref. 281, we will shortly comment on three 
points to illustrate how restrictive the FN mechanism actually is: 



• No LR symmetry: 



26,2013 3:10 WSPC/INSTRUCTION FILE keV'Merle 



40 Alexander Merle 

Left-right symmetry 90 ' 224 ' 289 ' 290 is by itself very restrictive: the RH charged 
leptons and neutrinos are members of one and the same doublet, which im- 
plies hi = gi, the symmetry would also force the LH and RH fcrmion dou- 
blets of each generation to have equal FN charges, /, = g t . Apparently, with 
such a strong restrictions, no semi-realistic pattern as the one in Eq. (50) can 
be obtained. This problem is particularly interesting since one mechanism 
to produce a suitable abundance of keV sterile neutrinos, namely thermal 
overproduction with subsequent entropy dilution, was mainly studied in 
LR-symmetric settings. 53,91 

• Difficulties with 5*0(10): 

The FN mechanism is typically applied in the context of Grand Unified The- 
ories (GUTs). In the case of SU(5), the RH neutrinos are total singlets, 
which leaves their FN charges perfectly unconstrained, just as required to 
achieve the desired RH neutrino mass patterns, cf. Eq. (51). However, in 
theories based on 50(10), the RH neutrinos are members of a 16^ repre- 
sentation, for each generation i. 291 Hence, their values would be strongly 
constrained, which destroys the freedom needed to generate a strong hier- 
archy among the sterile neutrino masses. 

• No need to run: 

Sometimes, non-fitting FN-based models are claimed to have a better match 
with the data when the couplings arc evolved to low energies by renormal- 
ization group running, especially in cases where the models under study do 
not lead to a good fit with data. 283 While this is a natural thought, since 
the predictions obtained by the FN mechanism should only be strictly true 
at a high energy scale, the running of the couplings is nevertheless fully 
negligible for low scale seesaw mechanisms as required by typical scenarios 
for keV sterile neutrinos. This conclusion could only be altered by having 
strong degeneracies between two or more light neutrino masses, at a relative 
level of roughly 10 -5 . However, since the FN mechanism generically leads 
to hierarchies and not to degeneracies, such effects can hardly appear in 
the models under consideration. 

As we have seen, the FN mechanism is very well suited to obtain scale suppres- 
sions and to generate hierarchies between masses. However, we have also seen that, 
although FN is actually more predictive than one might naively think, it is not 
powerful enough to predict complicated mixing patterns by itself, without varying 
coefficients of 0(1). An obvious way out is to combine the FN mechanism with an 
additional discrete flavor symmetry, which is the next point to be discussed. 

5.2.2. Mixed FN models 

We have now seen how one can build a model based exclusively on the FN mech- 
anism, apart from an auxiliary Zi symmetry that is "only" important if manifest 



26,2013 3:10 WSPC/INSTRUCTION FILE keV'Merle 



keV Neutrino Model Building 41 

CP violation is desired. The alternative type of models on the market uses the FN 
mechanism only as one among several ingredients, typically in combination with 
a non-Abelian (i.e., non-commuting) discrete flavor symmetry This approach has 
the great advantage of being able to exploit the powerful FN suppressions, while 
at the same time making use of the structure of the discrete group in order to pre- 
dict more sophisticated mixing patterns in the lepton sector. However, as always in 
model building one pays a price, which is in the case at hand the loss of minimality. 

The models on the market make use of an A4 symmetry, along with an additional 
small auxiliary symmetry Z 3 ° in combination with the FN mechanism. 273 ' 279 The 
model we arc going to discuss in some detail here is the A4 seesaw model with 
one ke V sterile neutrino (in the version with the other two sterile neutrinos being 
considerably heavier, with masses of at least a few GeV) presented in Ref. 273. This 
model is maybe the prime example in the literature to illustrate how to combine 
the FN mechanism with a non-Abelian symmetry in a smart way. in order to obtain 
a model with a keV sterile neutrino. 

When looking into the literature, typically the first piece of information given 
about a model is the particle content as well as the charge assignments, although 
obtaining this knowledge typically requires a lot of work and failed attempts. For 
the model under consideration, this information is presented in Tab. 2. We first 
of all see that, besides the necessary lepton doublets and singlets, two Higgs fields 
H u ,d are needed in the model, as well as seven fiavons <p, ip', ip", £, and 0. 

In addition, Tab. 2 gives us the information of how the fields transform under the 
different symmetries, either by directly indicating the representation [for SU(2)l 
and for A4] or by indicating the charges [for Z3 and E/(1)fn]- Note that the field 
is a total singlet under all symmetries but /7(1)fn, so is nothing else than the 
FN flavon. The other fiavons, however, are triplet or singlet fiavons, under A4, and 
they induce most of the structure. Finally, the sole purpose of the Z 3 assignment 
is to switch off certain terms that would not be desired, for example because they 
could spoil the mixing pattern. 

Next, we write down the leading order Lagrangian, where only terms that are 
singlets under all symmetries are allowed. Glancing at Tab. 2, making use of the A4 
group theory (cf. Sec. 3.2.3), and realizing that ui 3 = 1, the most general leptonic 
Yukawa coupling Lagrangian is given by: 

£y = - y -^\*eH{ipH d L) x - ^\M(<pH d L) v - V -^^H d L) v , (54) 
~\ gi JhR{yH u L) x - j-\ 92 N^(ip'H u L) v > - ^X>*N^{tp" H u L)x 

where the notation (...)i rre p refers to the corresponding combination under A±. Note 



° Z3 = {e,a, a 2 } is the cyclic group of three elements. Hereby, e denotes the neutral element and 
a is a non-trivial element with a? = e. 
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that the FN mechanism has already been applied in Eq. (54), which means that 
the invariance under J7(1)fn is only visible once the different powers of the small 
parameters A = (0) /A are taken into account. Note further that the notation has 
been slightly changed in order to be consistent with the rest of this review, which is 
why the Lagrangian presented here looks a bit different from the one in Rcf. 273. p 

Table 2. Particle content and representation/charge assignments (in our notation) of the 
SU{2) L x A 4 x Z-i x C/(1) PN model with one keV sterile neutrino. 273 SU(2) L is the usual SM gauge 
symmetry, and Z3 is an additional auxiliary symmetry in order for CP violation to survive, similar 
to the Z2 used in Sec. 5.2.1. For SU{2)j, and for A 4 the representations are given, while for Z3 and 
t r (l)pN the charge is indicated. Note that ui = e 2,r *' 3 . 
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Again, one has to choose a certain VEV alignment. In the case at hand, this 
alignment can, e.g., be achieved in a supersymmetric context by adding further 
auxiliary fields, 261 but also other origins might be plausible. As we had mentioned 
in Sec. 3.2.3, we will not embark this topic in this review, as achieving VEV align- 
ments is nearly a discipline by itself. In some sense, this belongs to the "dirty" part 
of model building, since it typically involves an essentially untestable sector and 
the involvement of many fields that are typically shifted towards higher energies. 
However, an alternative view would be that there are indeed many possibilities to 
achieve a certain VEV alignment, and the mere number of working possibilities is 
indeed a bonus of flavor models. 

Ending this discussion to enter the actual physics, the alignment chosen in 
Ref. 273 for the first triplet fiavon field is (cp) = (v v , 0, 0), which leads to a charged 
lepton mass matrix given by 

fy e X 3 \ 
M e = ^ y,X , (55) 
\ yj 

where Vd = {Hd) is one of the two Higgs VEVs. Note the FN suppression, which is 
clearly visible in Eq. (55): although different assignments have been used, the actual 
pattern of the charged lepton masses clearly resembles the one obtained from pure 
FN models, cf. Eq. (53). 

Up to now, the FN charges 31,2,3 of the RH neutrinos have not yet been fixed. 
With a VEV u = (£), the mass of the lightest sterile neutrino is given by 

M 1 =w 1 u\ 2gi , (56) 



p Note also the typo in Ref. 273 in the third charged lepton Yukawa coupling. 
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where w% can always be chosen to be real and positive by absorbing any phase 
into the hclds. Taking g\ large enough rapidly decreases this mass down to the 
keV scale, while the parameter w\ is still left for finer adjustments. For example, 
in Rcf. 273 this charge is suggested to be g\ = 9, which leads to M\ « 1 keV for 
A w 0.1 and u ~ 10 12 GeV. However, one has to remark that such a large gi can 
also be problematic, as the FN charge actually enhances the active-sterile mixing, 
6 e i ~ ^2 A~ gl . Hence, one has to be extremely careful that a too large value of 
gi does not spoil the validity of the model. 

Furthermore, Ref. 273 states that this lightest sterile neutrino decouples from 
the seesaw formula due to the strong X-ray bound, 6\ > 10 -8 for Mi « 10 keV, cf. 
Sec. 4.1. While this is certainly true in some approximation, one can argue that one 
does not have to care about decoupling Ni, since the FN charges 51,2,3 in any case 
cancel in the seesaw formula, cf. Sec. 5.2.1 and Rcf. 281. In particular, the quantities 
u>i and u from Eq. (56) will actually enter the seesaw formula and cannot necessarily 
be neglected. Alternatively, one could argue that according to the seesaw practicality 
theorem, cf. Sec. 4.2, we can in any case make use of the seesaw formula. 

In addition, in the approximation of a fully decoupled Ni , one of the light neutri- 
nos must actually be massless at tree-level: the seesaw fair play rule, cf. Sec. 4.2, 
states that two massive right-handed neutrinos can lead to at most two massive 
light neutrinos which might contradict future experiments. Hence, even though the 
decoupling of Ni is true to some extent, it is not always a good approximation. 

Let us nevertheless follow the path of Ref. 273 and see where it is leading us. 
The remaining task is to determine the mass and mixing pattern. As often the 
case for FN-inspired models, normal or inverted neutrino mass ordering are both 
possible. Ref. 273 suggests to take the alignment (if') = (vL,v',v') for both cases, 
while (if") = (0,v'',—v'') and (ip") = {2v'^,—v'^,—v'^) lead to normal and inverted 
ordering, respectively. With these choices, one can straightforwardly write down the 
Dirac mass matrices, 




(57) 



where v 



(H u ), and the RH Majorana mass matrix, 




(58) 



Using Eqs. (57) and (58), one obtains the full 5x5 neutrino mass matrices for both 
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orderings. These matrices can be fully diagonalized to yield the eigenvalues 
0. 



NO 



NO 
1 2 

NO 
l 3 

NO 

M 

NO 



ml 0) fl - 3c?) 



'2 
(0) 



1-24) 



w 2 u'X 2F " 
w 3 u"\ 2F ° 



4 0) 



(1 - Set) 



24) 



1° 


= m<°> (1 - 


Vet), 




IO 
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= m { 2 ° ] (1 - 






IO 
3 








IO 
4 


= w 2 u'\ 2F - 


- mf ] ( 


1-3,?) 


IO 


= w 3 u"\ 2F * 


(0) 

— m\ 


(1 - 6ei) 



(59) 



for normal and inverted ordering, respectively, where we have used the abbreviations 



,(°) 



6y 2 v" 2 v? 



(0) _ 3y 2 v' 2 v 2 



m 



(0) _ 2y 2 v" 2 v 2 



w 3 u"A 2 ' w 2 u'A 2 
_ mv^vu _ yav^vu 

— i £ 2 = TTT^ ■ 

w 2 u A !»3S A 



W3U 



"A 2 ' 



(60) 



Again, one has to be careful in the choice of the FN charges 172,3 in order not to 
spoil the smallness of e\ j2 - Still, the above patterns are somewhat compatible with 
the experimental data, and they can be used to calculate predictions for neutrino 
observables. 

The corresponding diagonalization matrices are given by 



[/(NOJO) 



/A 73 \ / 

L _i L 

L J_ J_ n 

V6 V3 V2 

1 

V 001/ \o, 








-%/6e 2 
















-V2e 2 , 



ei 0,2e 2 \ 
ei -e 2 
ei £2,-62 





-0(e 2 ). 



(61) 



Note that the upper left 3 x 3-block of the leading contribution is tri-bimaximal , 



2(14 



C^TBM 



/ 2 1 n 

1 1_ 

" V6 V3 V2 

\ L J_ J_ 

\ V6 V3 V2 



(62) 



This yields #13 = in particular, which is now known to be excluded, cf. Eqs. (4) 
and (5). Although the leading order matrix ui N0 ' 10 ^ is corrected by terms of 0(ei) 
and higher, the first (and actually also the second) order contribution do not do the 
job of generating a large enough #13. Even higher orders are likely to be too small 
to yield a reasonably sized #13. However, it is argued in Ref. 273 that this problem 
can be overcome by including higher-order terms in the Lagrangian, Eq. (54). Since 
such corrections are unavoidable, this is a good way of solving the problem, and the 
reference further shows that relatively large values of 613 are in fact possible. 

Summing up, the model presented in Ref. 273 comprises a fully working setting, 
as long as the parameters are carefully chosen. Nevertheless, although FN-inspired 
models are seemingly perfect to explain keV neutrinos, choosing a too large FN 
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charge can lead to problems with the strong observational bound on the active- 
sterile mixing. 

5.3. Models based on flavor symmetries 

Apart from the FN mechanism, which is already a simple kind of flavor symmetry, 
one can also use other symmetries. In the literature, one currently finds models for 
keV sterile neutrinos based on two relatively simple symmetries, L e — — L T and 
Qe, although other possibilities are not excluded. The important property that the 
existing models based on flavor symmetries have in common is that the keV sterile 
neutrino mass is completely forbidden at leading order. Hence, the "natural" value 
of mass of the DM neutrino N\ is exactly zero, but higher order corrections lift it to 
a non-zero value, which clearly follows a bottom-up scheme, cf. left panel of Fig. 4. 
Furthermore, if the other sterile neutrinos A2.3 obtain a mass already at leading 
order, these types of models simultaneously yield an explanation for Mi <C M?.^,. 
Of course, this does not necessarily have to lead to the absolute scale of keV, but 
the models indeed do explain a certain structure (or hierarchy) in the RH neutrino 
sector. 

5.3.1. L e — — L T symmetry 

The most straightforward symmetry to use is, similar to FN, a particular U(l) sym- 
metry which essentially plays the role of a generalized lepton number. The symmetry 
under consideration is called L e — — L T , for reasons that will become clear soon. 
For simplicity, we will abbreviate J- = L e — — L T . Based on considerations in 
Ref. 292, this T symmetry was proposed in Ref. 293 (two RH neutrinos) and in 
Refs. 294, 295 (three RH neutrinos). In general, the symmetry leads to a spectrum 
(0, m, m) for light neutrinos in the limit of T being conserved. However, while for 
Ref. 293 this could be argued to originate from the seesaw fair play rule, in partic- 
ular Ref. 295 showed that the symmetry intrinsically predicts an analogous pattern 
(0, M, M) in the heavy neutrino sector where no seesaw is at work, if an additional 
fi-T symmetry is present. A first application of this symmetry (and of similar pat- 
terns) to the case of keV sterile neutrinos was presented in Ref. 276, and a recent 
more detailed study was performed in Ref. 277. 

The charge assignments of this model are illustrated in Tab. 3. Note that T 
is a special type of f (1) symmetry, which obviously has its name L e — — L T 
from the charge assignments. At first sight, this model looks considerably more 
economical than, e.g., the one in Sec. 5.2.2. However, the reason for this is that the 
exact breaking of the symmetry was not specified in Ref. 277. This may be seen 
as incompleteness or inconsistency. On the other hand, the effect on the spectrum 
discussed in the reference is general enough that it is applicable for more than one 
type of breaking. 

With these assignments, one can immediately write down all allowed terms in 
the Lagrangian, in the most general setting of seesaw type II, cf. Eq. (A. 9). However, 
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Table 3. Particle content and representation/charge assignments (in our 
notation) of the SU(2) 2 X T model. 277 SU(2) L is the usual SM gauge 
symmetry. For SU{2)i, the representations are given, while for T the charge 
is indicated. 
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we will for illustration focus on the RH terms first, which are given by 

£m = 



-M R z {N 1R yN 2R - M R A (N 1R y N 3R 

MgM R - 





\(N R Y 



M R 2 



N. 



it 



^ M R 



h.c. = -^(N R )-M R N R + h.c.(63) 



The key point is that this matrix must have (at least) one zero eigenvalue, as one 
can easily see by computing the determinant: 



det M R = O-O-O+Mg-O-Mg+Mg-Mtf-Q-Mtf-O-Mtf-O-O-O-O-MX-Mg = 0. 

(64) 

Indeed, the eigenvalues of the matrix M R turn out to be (0, +M, —M) so that one 
sterile neutrino is exactly massless, while the other two are degenerate with masses 
M = V(M^) 2 + (M«) 2 . 

Applying the same logic to the full 6x6 neutrino mass matrix, one obtains 



C v = - 



h.c, where M v = 



( 



eii 



in, 



in 





r3 

m D rri jj 



M R 2 

M R 









\ 


u2 
m D 


u.3 
m D 




ml 2 


™ t d 




M R 


M R 3 


















/ 



(65) 



and \I> = (v\l, ^3L, N£ R , N^ R , N^ R ). In the seesaw type II limit, rn^f 
■m°£ <C M R , this matrix can be analytically diagonalized under the assumption 
mfj ~ mo- Even for this case one obtains two zero eigenvalues, one in the LH 
and one in the RH sector, while the other two pairs of eigenvalues are given by 

A± = ± yj { l^~ A { lf^ m D + 0(m 3 D ) and A± = ±V(M^) 2 + (Mg)* + 0(m%) for 

the light and heavy neutrinos, respectively. We can see that we furthermore need 
M R 2 ~ M R 3 to a good precision for the "small" eigenvalues X± to be small enough. 
This, together with the assumption on the entries of mrj, is in effect nothing else 
than introducing an additional fi-r symmetry, 296 in accordance with Ref. 295. 

While the eigenvalue spectrum looks promising, in the sense that a strong hi- 
erarchy in the heavy sector is generated, we still need to justify a non-zero sterile 
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neutrino mass M^.^ At the same time, the diagonalization of the light neutrino mass 
matrix is achieved with the so-called bimaximal mixing matrix, 



Ubm = -I I 775 . (66) 




which predicts the non-acceptable values d\ 2 = § an d 6*13 = 0. This is a general 
feature of T symmetry. 297 Both these problems can be cured by breaking the sym- 
metry. The exact mechanism of this breaking has not been discussed in Ref. 277, but 
instead several soft breaking terms, i.e., terms which suitably break the symmetry 
but do not introduce dangerous quadratic divergences, have been assumed. How- 
ever, the point is that no matter which breaking pattern is assumed, the resulting 
effect is very generic. 

Before we show how to cure the mass and mixing pattern by this breaking, let 
us note that by the Goldstone theorem 242,243 breaking of a global and continuous 
symmetry would lead to a massless Goldstone boson, r which is phenomenologically 
problematic. Possible ways out would be to gauge the symmetry, which would gen- 
erate a new massive gauge boson, or to use the discretized version of the symmetry 
instead: since only the U(l) charges ±1 are used in Tab. 3, one could instead use, 
e.g., a cyclic Z4 = {±l,±i} symmetry. The reader is invited to verify that the 
change of assignments from (1, —1, —l)u{i) to (i, —i, — £)z 4 f° r the first, second, and 
third generation leptons exactly reproduces the mass matrix given in Eq. (65). 

We will now turn to a solution of the problem of massless neutrinos. In Ref. 277 
it was suggested to introduce soft breaking terms s ^ e ' wt ' rr m^f and g e ^i L i L - TT ^ 
on the diagonal part of the matrix, 

M v -+ M v + diag( s | e , sl\ sY, S%, Stf, S R T ). (67) 

These new terms must be smaller than the "old" ones, since they arise from sym- 
metry breaking. In fact, even when viewed as explicit breaking terms, the whole 
consideration of a symmetry makes only sense in this limit. s The effect of these 
breaking terms can be studied analytically in the limit s£ Q ~ s, ~ S, and 
M R 2 ~ M R 3 ~ Mr, which allows to compute the eigenvalues of the light neu- 
trino mass matrix m„ = mj, — moM^m^ to be {A' + , A'_, A s , A' + , A'_, A s }, where 

A' ± = s±V2 \m L - =|] + r ^ M f+G (|0 , X s = s, A' ± = S±V2M R , and A s = S. 
Indeed, S < M R justifies taking S = O(keV), while M R can be of O(GeV), or 
larger. In fact, using the experimental data on the mass square differences, one can 
use these results to fully predict the light neutrino mass spectrum to be m\ = s + b, 



q In this model, it is actually the third generation neutrino which is massless. Hence, for the light 
neutrinos an inverted hierarchy will be predicted, while in the heavy sector it is A^3 which plays 
the role of the keV sterile neutrino. 

r In the case at hand, the boson would actually be a so-called Majoron. 2 ®^ 300 

s This is very similar to the isospin symmetry between protons and neutrons in nuclear physics. 
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m-x = s — 6, and m 3 = s, where b = 0.0489 eV and s = —3.9 x 10 4 eV, which yields 
\mi\ = 0.0486 cV, |m 2 | = 0.0494 eV, and |m 3 | = 0.0004 cV. This, in particular, 
confirms that inverted ordering is predicted. 

It remains to show that also the neutrino mixings turn out to be okay. Sticking 
to our simple example, it has been shown in Rcf. 277 that indeed the diagonalization 
matrix for light neutrinos can be modified compared to Eq. (66). First of all, the 
light neutrino mass matrix is given by 



/2c 22 c2 2 2 c2 , 

. rn n b m D m D b m D m D b s 

S + 2Mj mL ~ T1r~ 2M|~ mL ~ TTr~ 2M|~ 
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,S 1 m n S 



\" lL Mr 2M| "T7g~ * T TTf , 

This modifies the neutrino mixing matrix from bimaximal to 
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where e = : / m - \ ■ ( 69 ) 



16M 



However, this matrix still has to be combined with the mixing coming from the 
charged leptons. A suitable charged lepton mixing matrix is given by 301 

'l-A 2 /2 A A 3N 
Ui= \ -A 1 - A 2 /2 A 2 ) + C(A 4 ), (70) 
A 3 -A 2 1 . 

where A ~ 0.20 is the deviation of #12 from ^. Building the PMNS matrix C/pmns = 
U}U V in this basis (with non-diagonal charged lepton mass matrix), one obtains for 
the mixing angles 

(tan 2 12 , \U e3 \, sin 2 29 23 ) ~ (l - 2^2A - 2^2A 3 + 4A 4 , -j=, 1 - 4A 4 ^) , (71) 

or (sin 2 6*12, sin 2 6»i 3 , sin 2 6*23) =i (0.295,0.020,0.460), which is all consistent with 
data at the 3a level, cf. Eq. (5). 

Finally, it remains to be shown that the charged lepton squared mass matrix 
necessary for Eq. (70) to arise, 

/m 2 + m 2 A 2 m 2 A \ 
M,M/~ m 2 A m 2 +G(A 3 ), (72) 
V m 2 / 

can be obtained from softly broken T symmetry. Indeed, if the symmetry preserving 
charged lepton mass matrix as derived from the assignments in Tab. 3 is extended 
by soft breaking terms, s"' 3 , which would be zero in the symmetry preserving limit, 
one obtains 

b i b ' 





I mf e 








Mi = 





LLLL 

m i 


LIT 

m i 




V 


TU 

m, 





mf r . (73) 
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This yields Eq. (72) for the identification s e ^ = mf^X, s\ r = 0, m\ e = m e , mf^ = 
m M , and mj T = m T . 

The corresponding mass shifting scheme is depicted in Fig. 6, and it clearly re- 
sembles the bottom-up scheme from the left panel of Fig. 4. As illustrated in the 
figure, the masses of the two heavier neutrinos are exactly degenerate in the symme- 
try limit, and the lightest sterile neutrino is exactly massless. Once the symmetry 
is broken, the additional small soft terms will not only lift the degeneracy between 
the two heavier states, but they will also give a small mass, of the order of the 
breaking terms, to the lightest state. Hence, this scenario motivates a mass pattern 
Mi = O(keV) < M a , 3 . 

5.3.2. Qq symmetry 

Another proposal in the literature is a model for keV sterile neutrinos based on the 
symmetry group Qe, which was presented in Ref. 302. Q§ is nothing else than the 
double cover of the dihedral group D 3 = S^,* as explained, e.g., in Refs. 304, 305, 
306. 

The necessary group theory of Qq is outlined in Ref. 302, but we will comment 
on a couple of things in passing. This group is not much more difficult than A4: 
it has twelve elements, and four different one-dimensional irreps 1, 1', 1", 1"', 
as well as two different two-dimensional irrpes 2, 2'. Note that 2 is a pseudoreal 



t D% is the full symmetry group of a triangle, i.e., three rotations and three reflections, which 
amounts to a total of six elements. This is why some authors use the notation "Dg" for this group, 
which is unfortunately a bit ambiguous as it could also be interpreted as the full symmetry group 
of a hexagon. 303 

I M 2 =M 3 ;>GeV 
▼M 2 £GeV 

L e -L^-L T & /i-r 



r Mi ~ke V 

* M,=0 

Fig. 6. Mass shifting scheme for softly broken T symmetry. Breaking the symmetry lifts the 
sterile neutrino mass M\ to a non-zero value, which clearly resembles the bottom-up scheme, cf. 
left panel of Fig. 4. (Figure similar to Fig. f .a in Ref. 281 ) 
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representation, i.e., it contains only matrices with real entries, but the squares of 
some of them are nevertheless negative. An important difference to At is, however, 
that Qq has no three-dimensional irreps, so that a model based on Qq can unify at 
most two of the three known generations of matter. Because of this fact, one might 
criticize the aesthetic power of this group, but one could also argue that part of the 
data actually seems to favor a unification of only two generations of leptons. For 
example, when looking at charged lepton masses, and m T are much closer to 
each other than they are to m e . One can always take on a very pragmatic point of 
view and judge a model by its consistency with experimental data rather than by 
requiring it to predict certain patterns that may or may not exist in Nature. 

The model under consideration is relatively economic since only three singlet 
flavons S x , y ,z and one doublet flavon D arc introduced in addition to the lepton 
sector. The detailed charge and representation assignment is outlined in Tab. 4. 
Also in this model, two auxiliary symmetries have been introduced in order to 
suppress certain terms in the Lagrangian. 

Table 4. Particle content and representation/charge assignments (in our notation) of the 
SU(2)l x Qe x Z>2 X %3 model with one keV sterile neutrino. 273 SU (2) l is the usual SM gauge symmetry 
and Z2,'i are additional auxiliary symmetries. For SU(2)r, and for Qq the representations are given, 
while for ^2,3 the charge is indicated. Note that u> = e 27 ™/ 3 and that 2* = 2 in Qg. 
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Using the representation content as well as the group theory of Qe, we can again 
write down the lowest order Lagrangian for the charged leptons," 

Cl = -^(l^HLD^Sl-^i^HLD^Sl-^ie^HL^riSlf + h.c. (74) 

Here, we have made use of the composition rule for two doublets 2 = (xi, X2) T and 
2' = (yi,V2) T , whose tensor product 2® 2' can be decomposed into two singlets, 
1" = {xiiji - X2JJ2) and 1"' = (xxyi + .T21/2), and one doublet 2 = (x 2 yi, xiy 2 ) T ■ 
Furthermore, we have used 1' (g) 1' = 1 and 1" (g) 1" = 1"' <g> 1"' = 1'. Note that the 
whole Lagrangian in Eq. (74) arises as correction to the leading order, and hence 
all the terms inside are non-renormalizable and suppressed by a high mass scale A. 
Assuming singlet flavon VEVs (Si) = Si for i = x,y, z, the charged lepton masses 



"Note that this is the Hermitean conjugate of the Lagrangian presented in Ref. 302. Furthermore, 
we correctly take the mass terms to be negative, as is conventional for a Dirac mass, although one 
could simply redefine the couplings to account for that. Finally, since we are effectively taking the 
complex conjugate, one has to keep in mind that 1"* = 1"' in Qq. 
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obtained after clcctrowcak symmetry breaking are given by 

/ s x \ 2 1 1 

m e = {^-j-J Y e v, m,, = -r(Y x s x + Y y Sy)v, m T = —(Y x s x -Y y s y )v. (75) 

Indeed, a first tendency that had been anticipated reveals itself in these eigenvalues: 
the representation content was chosen such that the muon and the tau unify in 
doublets of cf. Tab. 4, while the electron transforms as a singlet. Accordingly, 
the masses of \i and r in Eq. (75) are identical up to a splitting 2Y y s y v, while the 
electron mass is suppressed (compared to m^. T ) by another power of the high energy 
scale A. 

The corresponding Lagrangian for the neutrino sector turns out to be 

C» = —(NElf^HL! - ^-(N^HL D ) 2 ,D* - ^ (N~dHL* d )±S x 

—^(NDi^HL^S*) 3 - -^(N^HL D ) v S*S y S* z (76 

~[W^F^d]i - I^WdFNdW [{D*)\, \^mT c NAvS* x S v + h.c 

where the RH Majorana masses are taken to be real. Note also the missing factors i 
in the corresponding equation in Ref. 302. In Eq. (76), we have made use of further 
tensor products, such as 1" ® 1"' = 1 or 2' ® 2' = 1 © 1' © 2'. The reader is invited 
to work out the exact form and the invariance of the remaining terms in Eq. (76) 
her or himself, using the information on the group Qq given in Ref. 302. 

The most important aspect of Eq. (76) can be found in the last line of the equa- 
tion: since No transforms as a doublet 2' under Qq, the combination (Nd) c Nd 
indeed contains a trivial singlet 1. To be precise, the tensor product of two doublets 
2' = (xi, x-i) T and 2' = (yi, y2) T contains a term 1 = (xiU2 — x^yi), and the conju- 
gations inside (iVi) c do not matter since 2' is real. In conclusion, the corresponding 
term in the Lagrangian yields a tree- level mass for N% 3, 

- l(N^fN D ] 1 + h.c. = -^- (7V2iV 3 c - AyVI) + h.c. (77) 

However, the first RH neutrino, N\, transforms as a non-trivial singlet and, similar 
to what happens to the electron in Eq. (75), its mass is suppressed by two more 
powers of the high scale A compared to the masses of ^2 3. A natural hierarchy 
Mi <C A/2, 3 in the sterile neutrino sector is predicted, and hence we have a good 
reason for taking Ni to be the keV sterile neutrino. 

Taking the VEV of the doublet fiavontobe (D) = [d\, c?2) T , one can immediately 
read off M R from Eq. (76), 

/0 

M R = M a |+-| M h d{ 
\0M a 
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whose eigenvalues are given by Mi = ^j^-M c and M 2i3 ~ M a =p $s(df + <$£). The 
Dirac mass matrix, in turn, is given by 

/0 ad 2 adA {8sl \ 

= V I fidi ^s z J « + -p es x s y s x u. (79) 
\0 7s z /3d 2 / \ -es x s y s z / 

Due to Mi <C M 2i 3, the authors of Rcf. 302 chose to integrate out 7V 2j 3 only, although 
they could actually have applied the seesaw formula according to the keV seesaw 
practicality theorem, cf. Sec. 4.2. Restraining from that, according to the seesaw 
fair play rule, cf. Sec. 4.2, one light neutrino mass eigenvalue will be precisely zero. 

Applying this procedure, the authors arrive at a 4 x 4 neutrino mass matrix 
given by 

»r - ("p . (so) 

where A = (S'v, 0, 0) T , and the upper left 3x3 block is 
, / 2a' 2 «'(/?' + Y) «'(/?' + 7') 



M 3X3 = J_ 
M a 



a'(P'+i) 2/3'i /3' 2 + 7 ' 2 
\a'(p'+i) f3 ,2 +i 2 2/3V 

a'(2/3'-j') -a'(2/3'- 7 ') N 



+e d — a>{2p> - i) 2p'i | + O(cj), (81) 

1 -a'{2p'-i) -2p'i 

where d = dl + d2 and e d = rfl 2rf d2 , such that = in the limit d\ = rf 2 . Furthermore, 
the abbreviations used are a' = aj, /?' = /3j, 7' = 7^, and (5' = S (x) 3 - Note 
that the latter definition suggests that the parameter 5' is actually quite small. 

In Rcf. 302, the 4x4 neutrino mass matrix is approximately diagonalized in the 
limit of small 6', which is done by a matrix £/4 X 4 given by 

cos(9e v sin6»e v 

i(l-\R®W) R \ (Va 0\ , ^ I sine w co S « w 1 



\/2 V2 

(82) 

and i? = (5'v/Mi, 0, 0). Furthermore, the abbreviations p = arg(a'), er = arg(/3' + 
7'), and sin 9 = . ^ a ! have been used. Note that Va is exactly unitary and 

Vl/3'+7'l 2 +2|a'| 2 

[^4x4 is approximately unitary for small \R\, which indeed requires a relatively tiny 
5' [in Rcf. 302 it is taken to be of 0(1O -13 ' 5 )]. Note, however, that the expression 
given for Va only diagonalizes the leading part of M 3x3 , i.e., the one for which 
e d = 0. 

The mass eigenvalues of the light neutrinos resulting from the diagonalization 
of M 3x3 from Eq. (81) are given by 

mi =0 , m 2 = — w , m 3 = — v , (83) 

iK£ a lVl a 
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where indeed one light neutrino is exactly massless, just as anticipated. In turn, 
this model predicts not only normal mass ordering, but in fact a normal hierarchy 
among the light neutrinos. 

The authors also calculate the neutrino mixing angles resulting from Eq. (82). 
They argue that the phases p and a are unphysical, which is correct, but still they 
appear in their expression for 623- This is incorrect, as one could have seen when 
following the detailed diagonalization procedure of the leptonic mass matrices as 
outlined in the Appendices B of Rcfs. 262, 307. However, one can easily cure the 
results for the mixing angles as obtained in Ref. 302 by simply setting p = a = 
in their epxressions. Then, since U e = 1 due to the charged lepton mass matrix 
being diagonal, the physical mixing angles are identical to the mixing angles from 
the neutrino sector alone, Oij = Q\a . Their lowest order expressions are given by 



sin ~ sin 9 , tan 623 



4/3'V 



sin t/13 



V2a'(2/3' - 7') 



(84) 

where corrections of 0(e^) have been neglected. Note that the angle 6*13 as well 
as the deviation of 6*23 from ^ arise only as higher order corrections. Nevertheless, 
we now know that the parameter of e<2 must be large enough to yield an angle $13 
as large as measured, cf. Eqs. (4) and (5). However, the mixing angle #12 is not 
predicted and must be matched to the data. 

Let us end this section by an interesting feature of the model, which is illustrative 
for how flavor models could be tested. Comparing the equations for 6*13 and 6*23 in 
Eq. (84), one can derive a correlation between these two angles: 

2v^2/^''y' 

tang 23 ^l- k( ^_ y)| S ingi 3 . (85) 

This equation implies that 623 = 5 if 6*13 = 0, but also that #23 7^ § if #13 7^ 0. 
More precisely, the sign of /3'^' determines whether 6*23 is smaller or larger than ? . 
This is particularly interesting when looking at the newest global fits, cf. Refs. 15, 
16, 17. This correlation has been numerically verified and matched to the actual 
data in Ref. 302. Such correlations are the keys to identify whether a certain model 
finally turns out to be close to reality, or not. 

After having discussed several models based on flavor symmetries, we will now 
turn to other ways to justify the keV scale for sterile neutrinos. A completely differ- 
ent mechanism to generate very small mass scales, based on the existence of extra 
spatial dimensions, will be presented in the next section. 



5.4. Models based on the split seesaw mechanism 

A somewhat orthogonal approach to explain the keV scale is by mechanisms that 
are not based on some flavor symmetry. One possibility is to have a non-trivial 
space-time geometry, and the most generic approach is making use of extra spatial 
dimensions. The corresponding mechanism which can explain the existence of a keV 
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scale is the so-called split seesaw mechanism. 308 Unfortunately, we will here not be 
able to review all the theory behind extra dimensions. Instead, we guide the reader 
to excellent introductions to and reviews of extra dimensions, such as Refs. 309, 
310, 311. However, the amount of knowledge on extra dimensions that is necessary 
to understand the key concepts of the split seesaw mechanism is relatively limited, 
and we will in addition give brief explanations of the required concepts in passing. 

5.4.1. The split seesaw mechanism 

The split seesaw mechanism has been proposed in Ref. 308. It is based on a 5- 
dimensional (5D) model compactified on S 1 jZ-i. the extra dimension is rolled to a 
circle S 1 , parametrized by an angle <j> £ [0, 27ri?]. In addition, the pairs of coordinates 
(cj>, 2ttR — </>) are identified with each other, which is called orbifolding. Hence, the 
5D coordinate can be taken to be ?/ £ [0,1], where I = ttR. The model contains 
two branes (= 4D subspaces), one at y = to which all SM fields are confined and 
another hidden brane at y = I. Such a setting leads to a reduced 4D Planck scale 
Mp, which is related to the fundamental Planck scale Mq by Mp = M$l. Note that 
this means Mp < Mq, due to l~ x < Mq. 

The next step is to write down the 5D action of the model for a 5D Dirac spinor 
= (^i, *&r) t with a bulk mass m in five dimensions: 

i 

S = J d*x J dy M (i^T A d A ^ - m*#). (86) 
o 

Here, the fundamental mass scale Mq serves the purpose of correctly normalizing the 
unit of S. The index A runs over 0, 1, 2, 3, 5, where 5 labels the extra dimensional 
coordinate y, and the corresponding gamma matrices are defined as 

r» = r , r5 = (J • ( 8? ) 

The important point is that the LH and RH projections are defined by the new T 5 
matrix, i(l T r 5 )$. v 

The next step is to integrate out the extra dimension, in order to obtain an 
effective model in 4D. A typical feature of such settings is the existence of a Kaluza- 
Klein (KK) tower, which is a stack of 4D fields with identical quantum numbers 
(up to chirality) but increasing masses. A suitable ansatz for a decomposition in 
KK-modes in our case is given by: 312 

*L, R (x»,y) =E<«(^) /£&(V). ( 88 ) 

n 

where the index n labels the different KK modes. Note the separation of the coordi- 
nates on the right-hand side of Eq. (88). This ansatz has to be inserted into Eq. (86), 



v See Ref. 310 for a discussion of the issue of defining chirality in theories with extra dimensions. 
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where we should be careful to note that Wc^* = * lu^B^ L + ^ R^d^ R while 
*"r 5 5 5 * = Wr^(-il)d 5 ^ R +W^(+il)d 5 '^L, since T 5 anti-commutes with all T M but 
of course commutes with itself. Inserting the KK-decomposition, one obtains 

s = j dAx J2 [v^V^f] (/ d y M o f { L ) \y)f i L k \y)) 

dyM f R n) *(y)f R k \y) 



+ / d+x ^Vi"W (fdy Mo f L n) \y){-d y - m )f R k \y) 

J n,k \ J 



J^x^rW (JdyMof^ 

n,k ^ 



(y)(+d y -m)f L k \y) 



(89) 



In order to yield the correct 4D interpretation, this action should be equivalent to 
S= fd A xY, (i¥^d^ { L n) - m n ^/4 n) + (L^R,W^*)), (90) 

J n 

which will be fulfilled if the following two relations hold, 



dyM f L n l*(y)f L k Uy) = S nk and (±d v - m)f L n) R (y) = m n & R (y). (91) 



For the zero mode, defined by too = 0, the bulk profile is obtained by solving the 
second Eq. (91), and it is given by a constant C times a function of y: 



f^ R {y) = Ce^ m ^ where C 



2to 



1 



,2m! 



IVMo 



(92) 



The correct normalization, cf. first Eq. (91), is fixed to yield a canonically normalized 
4D kinetic term. The solution given in Eq. (92) makes it right away evident that, 
for the RH solutions, the function value at the hidden brane (y = I) is much larger 

than the one at the SM-brane (y = 0), since , m — - = e 3> 1. This observation 

will be crucial in the what follows. 

In order to obtain the split seesaw mechanism, one applies the above procedure 
to a 5D RH-neutrino field whose zero mode is purely right-handed and which shares 
at the SM brane both, a Yukawa interaction Xi a = 0(1) and a baryon number minus 
lepton number violating Majorana mass k^vb—l •> 



S = d*x / dy 



M (*g*T*a A *g 



,(0) ; 



(93) 
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where i = 1, 2, 3 is a generation index. Applying Eq. (92), one obtains the effective 
4D sterile neutrino masses and Yukawa couplings, 



, r vb-l 2nii . A iQ / 2m l _ - / Mj 

Mi = Ki M n eW_i ' A » = ^V^^T = Aj ^ /rT:— - (94) 



KiVB-L 

Now we can see what happens: the sterile neutrino masses Mj receive two sup- 
pressions in Eq. (94), one from the ratio between the Majorana (B — L) scale and 
the fundamental Planck scale, VB M( ^ "C 1, and another even more significant one 
from the term with the exponential, which suppresses the natural mass scale rrn 
as e~ 2niil for Im.i 3> 1. It is this latter exponential suppression which makes the 
split seesaw mechanism so strong: a relatively mild mass splitting mi < 11J2 < 7713 
is enhanced by the exponential factor to yield physical masses Mi <C M 2 <C M 3 . 
The corresponding mass shifting scheme is depicted in Fig. 7. Clearly, this type of 
suppression follows the top-down scheme, cf. right panel of Fig. 4, just as FN type 
models do. 

An interesting feature of the split seesaw mechanism is that the seesaw is guar- 
anteed to work, similar to what happens in FN inspired models, cf. Sec. 5.2. Indeed, 
calculating the seesaw formula like in Eq. (13) with v = (H°), one obtains 

(muU = E A « A *^ = E ^Z-l (95) 

i i 

Hence, apart from a very minor enhancement proportional to k" 1 , the Majorana 
scale contribution oc Vb—l suppresses the electroweak scale contribution cx v. Note 
in particular that Eq. (95) allows for any light neutrino mass spectrum, may it be 
hierarchical or quasi-degenerate: the strong hierarchy in the sterile neutrino masses 
docs not translate into the active neutrino sector. 



\ Af 3 > Af 2 

Split \ 

Seesaw \ M 2 ~10 n GeV 



2 m; 



Mi~keV 



Fig. 7. Mass shifting scheme for split seesaw. The wave functions in the extra dimension suppress 
the physical sterile neutrino masses, which clearly resembles the top-down scheme, cf. right panel 
of Fig. 4. (Figure similar to Fig. l.b in Ref. 281 ) 
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Note further that there is less of an enhancement of the active-sterile mixing in 
the split seesaw mechanism than in FN inspired settings (cf. Sec. 5.2.2): 308 



fl ^3-10-' p t a S lal igl ■ (96) 




Still, a moderate enhancement is present: naturally one would expect 9\ ~ ^jP- oc 
Mf 1 , while fr om Eq. (96) we have the slightly milder relation d\ cx M x ' . 

The main advantage of the split seesaw mechanism is, as already mentioned, that 
a small hierarchy in m, translates into a large hierarchy in Mf. for example, taking 
mi ~ 242- 1 and m 2 ~ 2.3/- 1 leads to M x /M 2 ~ 10~ 17 , or M x = 0(1 keV) for 
A/2 = O(10 n GeV). However, one still has no explanation for where the hierarchy 
within the m, originates from. Furthermore, one could argue that a whole new 
hidden sector - featuring a lot of new particles from the KK-tower that can hardly 
be detected - is introduced to explain nothing more than the mass splitting of the 
sterile neutrinos. While the latter issue is a general problem in model building and 
there is hardly any way around it, the former problem can be circumvented by 
extending the split seesaw mechanism by a flavor symmetry. This is exactly what 
we will do next. 



5.4.2. The A4 extended split seesaw 

The up to now only extension of the split seesaw mechanism by a flavor symmetry 
was presented in Ref. 313, and it makes again use of our "good old friend" A4. The 
authors of this reference made use of the systematic classification of A4 assignments 
presented by Barry and Rodcjohann in Ref. 314. In this classification, ten different 
assignment schemes, called A, B, J have been identified. Out of these, some do 
not even contain RH neutrinos and are hence not very useful for the case at hand. 
Others assign the three RH neutrinos to a triplet 3 representation of A4, and also 
these schemes are not useful since they would lead to a degenerate mass pattern, 
M R {N R ) C N R = M R [(N 1R ) C N 1R + {N 2R ) C N 2R + {N 3R ) C N 3R ]. Another assignment is 
ruled out by the experimental values of the mass square differences, and out of the 
remaining two possibilities, the simpler one has been chosen for the discussion. 313 

The resulting model uses the same assignment (scheme E) as the original model 
from Ref. 315, but in addition to the Z 2 parity there is an additional Z5 auxil- 
iary symmetry imposed, which is important to constrain the fiavon potential. The 
resulting representation and charge assignments are detailed in Tab. 5, where in 
addition to the lepton sector and to the SM-like Higgs, one triplet fiavon tp v t and 
three singlet flavons ipn (i = 1, 2, 3) have been introduced. 

These assignments yield a non-diagonal RH neutrino mass matrix, so that the 
couplings Ki in Eq. (93) have to be promoted to matrix components k^ . Then, the 
Yukawa coupling matrices for the Majorana masses of the RH neutrinos and their 
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Table 5. Particle content and representation/charge assignments (in our nota- 
tion) of the SU(2)l X A4 X Z2 x Z5 extension of the split seesaw mechanism. 313 
SU(2)l is the usual SM gauge symmetry, and ^2,5 are additional auxiliary sym- 
metries. For SU(2)l and for A4 the representations are given, while for Z2,5 the 
charge is indicated. Note that ui = e 2 "" 1 / 5 . 



Field 


■^1,2,3 


ei,2,3R 




N 2 


N 3 


H 








¥>3,Z 


SU(2) L 


2 


1 


1 


1 


1 


2 


1 


1 


1 


1 


A 4 


3 


3 


1 


V 


1" 


1 


3 


1 


1" 


1' 


Z 2 


+ 




+ 


+ 


+ 


+ 


+ 








Zs 


UJ 3 


OJ 4 


1 


1 


1 


1 


u; 2 


u; 3 


CO* 





bulk masses are given by 

'a 0\ 

6 , m = diag(mi, 7712, 1713). (97) 
^0 6 0/ 

With a general triplet flavon VEV (</Vt) = w 2 ; W3) , one obtains the couplings 
Xia, at higher order 



r y^i ( fl(i^)oo + y 2 ^^(L^)oo" + ^JVaH^W ~> Afa^A,^, (98) 



A 

which leads to 



1 



'jtf \ /l 1 1\ /in 



A = - I yl I I 1 lj u? I I u 2 I . (99) 
Inserting the SM-likc Higgs VEV, one obtains the Dirac mass matrix, 



<y\ \ /l 1 1 \ lu x 
j/£ 1 cj cj 2 u 2 
y v J \1 u 2 lj J V u 3> 



ml = jdiag(Cx,C 2 ,C 3 ) 



where Ci = y e2 ,n"'_ 1 . The RH mass matrix, in turn, is given by 

laC\ \ 
M K = diag(d, C 2 , C 3 ) k diag(d, C 2 , C 3 ) = 6C 2 C 3 (101) 

V bC 2 C 3 / 

which yields the eigenvalues 

Mx = aCfv B - L , M 2 = M 3 = bC 2 C 3 v B - L - (102) 

Note the degeneracy of M 2 3 in Eq. (102), which may be crucial for baryogenesis. 29 
The light neutrino mass matrix can be immediately calculated using Eq. (13), 



m„ = m£ = -X 1 k~'X = \ if \ \ yxy \ \ if \ (103) 
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where x ; 



(Vl) 



2v2 b V3 and y : 



(vl) 



b 



As to be expected, the small factors 



Ci cancel out completely from the light neutrino mass matrix, Eq. (103). Assuming 
U\ = ru and U2 = W3 = u in {<J> v ,t), one obtains suitable eigenvalues for the light 
neutrino masses, which can be brought into agreement with the measured values 
of the mass square differences. 313 However, the problem is that this case leads to 
#13 = (and also to #23 = 7r /4), which is experimentally excluded, cf. Eqs. (4) 
and (5). Even worse, the active sterile mixing turns out to be too large by about 
three orders of magnitude, 313 which completely destroys the validity of the scenario. 

The way out of this dilemma is to extend the model further, by also including a 
seesaw type II contribution, cf. Eq. (18). That is achieved by adding a Higgs triplet 
field T, cf. Eq. (16), as well as two more flavon fields as illustrated in Tab. 6. This 
leads to an additional additive contribution S^ 1 to the action, given by 



in 



i A xdy5{y) ( 



yi — + y * — 



1 



LicxoT^L + h.c. 



(104) 



and the Higgs triplet VEV vt leads to an additional LH neutrino mass matrix 




(105) 



where s = y^^ 1 



given by 



u/A 



and m' 



Thc final neutrino mass matrix is 



n' v tm! v 
tm' v sm! v — iym,, 
,tm' v ym v 



tm' v 
ym u 
' v - 2ym vJ 



(106) 



which still leads to #13 = and #23 = 7r/4. The authors have recognized this 
problem, and they have argued that next-to-leading order terms could be able to 
cure it. However, no convincing discussion had been presented. 



Table 6. Particle content of 
the type II extension in addi- 
tion to Tab. 5. 



Field 


T 


fv,s 


<f>v,t 


SU(2) L 


3 


1 


1 


A 4 


1 


1 


3 


Z 2 


+ 


+ 


+ 


Z S 


^ 2 


uj 1 


1 



Putting aside this problem, however, the seesaw type II extension can provide 
realistic expressions for the mass square differences and for 6*12, which in the limit 



26,2013 3:10 WSPC/INSTRUCTION FILE keV'Merle 



60 Alexander Merle 



m,, look like: 



Am| x ~ (2s-y)\/8i 2 + y 2 m 2 , 

(2s -y- ^WTf)' 



Am 2 31 



ml, 



tan 



(y + ^ 2 + y 2 )' 
8*2 



(107) 
(108) 

(109) 



Note that, assuming the parameters to be complex in general, normal and inverted 
ordering are both possible in this model. Furthermore, this extension overcomes the 
problem with the X-ray bound: the LH mass matrix, Eq. (105), does not contribute 
to active-sterile mixing, cf. Appendix A. 2. Hence, one can choose the problematic 
elements in mo, Eq. (100), to be small enough while still keeping a realistic light 
neutrino spectrum alive by the type II contributions. 

We end this section by pointing out that the VEV alignment used in this model 
can actually be derived from orbifolding, as described in Rcf. 313. Furthermore, 
the reference discussed also various possibilities for leptogenesis within the model. 
These discussions, however, are beyond the scope of this review. 



5.5. Models based on the extended seesaw mechanism 

The next mechanism which we want to discuss here is the so-called extended seesaw 
mechanism. This mechanism was first proposed in Ref. 316 in order to solve the solar 
neutrino problem, and it was later on revived in Refs. 279, 317 for the purpose of 
explaining keV sterile neutrinos. While this mechanism goes beyond the "minimal" 
picture of having only three RH neutrinos, it nevertheless yields a relatively natural 
explanation for the existence of keV sterile neutrinos, at least when applied to 
concrete models. 



5.5.1. The extended seesaw mechanism 

Let us start by discussing the mechanism itself. The basic idea is to introduce, in 
addition to the three RH neutrinos Njh, another chiral singlet field Sr which is also 
right-handed. w The next step is to postulate the following Lagrangian, 

Ces = -vZm D N R - JSnfM^NR - ^{NrYMrNr + h.c. (110) 



w Note that we depart from the notation used in Refs. 279, 317, since in these references it is 
not so transparent that the singlet field called S must actually be right-handed. If it was not, its 
left-handed counterpart Sr, would exist, which could then be used to form an undesirable Dirac 
mass term for the singlet, Cs,vira,c = —Sl^/id^Sr + h.c. To make clear that such a term does 
not appear, we decided to change the notation, S — > Sr. 
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Up to this point, the Lagrangian presented in Eq. (110) is indeed nothing more 
than a postulate. The reason for this statement is that the Lagrangian is in- 
complete: gauge invariance would not forbid a Majorana mass term for Sr alone, 
£s,Majorana = — \ {Sr) c Mm,sSr + h.c, but this term does not appear in the above 
Lagrangian. However, the problem can be cured easily, by imagining that the term 
(Sn) c Mg Nr could arise from the VEV / = (4>s) of some singlet scalar field 4>s, 
USN(S R ) c (j)sNR -> {S R ) c (y S Nf)NR. If one uses a discrete symmetry, e.g. Z 4 , one 
can give a charge of ±i (=pi) to Sr {<fis), m which case the mixed mass term 
(Sr) c MsNr in Eq. (110) would be allowed, while a Majorana mass term for Sr 
alone would be forbidden, at least at tree-level. Of course, other solutions of this 
problem are possible. 

Let us press on and take Eq. (110) for granted. Then, Ms is a 1 x 3 matrix, since 
there is only one Sr versus three N^. Hence, one can rewrite Eq. (110) in terms 
of a 7 x 7 mass matrix Mj x7 , similar to the generic way of rewriting the seesaw 
matrix in Appendix A.l: 



1 



m D 




: ES = ~(VE, (Sr) c , (Nr)-) I ml Mr Mj | | S R | + h.c. (Ill) 

M s 

The important point of this Lagrangian is that there is essentially no constraint on 
the size of Ms- In particular, there is no need for Ms to be as large as the size one 
would naturally anticipate for Mr, since the different masses do not have the same 
origin (as explained above, Mr is a bare Majorana mass while Ms is part of a mixing 
term). On the other hand, there is also no reason to have Ms < Tod = 0(v), which 
would typically lead to problems with neutrino oscillations. 317 Hence, the most 
attractive possibility is to assume the following mass hierarchy: 

m D < M s < M R . (112) 

In this case, we will see that the inverse seesaw can indeed motivate a keV scale. 

The first step is to integrate out the RH neutrinos from the Lagrangian in 
Eq. (110). In theory, this means that the RH neutrino fields are assumed to be 
very heavy, and hence not dynamical at low energies. Then, one can neglect their 
kinetic terms and solve their Euler-Lagrange equations explicitly, just as done for 
the heavy fermions in the FN mechanism, cf. Sec. 3.2.2. In practice, this amounts 
to a seesaw-like formula, and the remaining "low" energy 4x4 matrix is given by 

M^=f m f%f D T m ° M \ M l). (113) 
iMs^ 1 ) m T D MsM^MT) 



The determinant of this matrix can be easily computed by the formula det 



det (A) det {D — CA^ 1 B), which holds for block matrices with an invertible subma- 
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trix A. With (M^ 1 ) = M^ 1 , one obtains: 
det (A/ 4x4 ) = det (m^M^m^) 

xdct (MsM^Mg - M s M^mF D (m^)" 1 M R rn^m D M^ 1 Mf \ = 0. (114) 

Since the determinant of a matrix is the product of its eigenvalues, this means that 
this minimal version of the extended seesaw mechanism intrinsically predicts one 
light neutrino to be exactly massless. 

In the case of Ms S> mj, one can now apply the seesaw type II formula, cf. 
Eq. (A.16), to Eq. (113): 

m v = M 3x3 = moM^Mg (MsM^M^)' 1 M s M R l m T D - m D M^m T D . (115) 

Note that the right-hand side of Eq. (115) does not vanish. 317 The reason is that 

is a number instead of a matrix, due to the structure of 
the matrix Ms- Hence, one cannot rewrite it as (AfJ ) 1 MrM^ 1 . In short, Ms as 
a 1 x 3 matrix is not invertible. 

The above block-diagonalization of Af* x4 also yields a singlet mass eigenvalue 
m s given by 

m s = M s M^Mj. (116) 

This equation reveals the whole trick of the extended seesaw mechanism: although 
the matrix product in Eq. (116) looks like a seesaw type I formula, i.e., m s <C 
|Ms,i|, we had chosen Ms 3> mo which enforces m s ^> \(m v )ij\. Due to these 
hierarchies, we have obtained a motivation for m s = O(keV) while m v < O(eV). In 
particular, due to the suppression of the new scale Ms in Eq. (116), the extended 
seesaw mechanism is another example of the top-down mass shifting scheme, cf. 
right panel of Fig. 4. However, the difference here is that a completely new mass 
scale is suppressed, and one actually has four sterile neutrinos instead of three. 

Although the extended seesaw mechanism in principle yields an explanation for 
the appearance of a scale that could be potentially of O(keV), we have already 
remarked after Eq. (110) that the required Lagrangian is a bit unmotivated in 
general. A suitable way out is to enforce the structure given in Eq. (110) by a 
flavor symmetry and, not too surprisingly, the simplest example that is known in 
the literature is based on an A4 symmetry. 317 

5.5.2. The A4 extended extended seesaw 

The A4 extension of the extended seesaw mechanism 317 can resolve some of the 
problems mentioned in the previous section. Again the task is to find suitable charge 
assignments, which arc summarized in Tab. 7. Note that this table also contains the 
flavon fields 95, <p' , tp" , £, and x-i while Sr is the additional singlet fcrmion. Using 
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these alignments, the leading order Lagrangian of the lcptonic sector is given by, x 

+ l^(N m yN R1 + ^\2£,'{N m yN m + ^(NrsYNrs + h.c, (117) 

where A denotes the cutoff scale. Note that we have again reported the Hermitean 
conjugate part of the Lagrangian presented in Ref. 317, in order to conform with 
our notation. However, except for Xi and p, all couplings and VEVs can be taken 
to be real by absorbing the phases into the fermion fields, and some phases will be 
forced to be equal by the VEV alignment. By radiative breaking, one can obtain 
the following alignment: 318 

(<P) = (v<p,0,0) (<p') = (v v ,v,p,v v ), {f") = (0,-v v ,v v ), (£) = {£') = v v , (x) = u. 

(118) 

This leads immediately to a diagonal charged lepton mass matrix, M e = 
^ £ -diag(i/ e , yp, y T ), which means that v v yi must be real for i = e,fi,T. Note that 
the hierarchy in the charged lepton sector is not explained by this model. 



Table 7. Particle content and representation/charge assignments (in our notation) of the 
SU(2)l x A4 X Z4 extension of the extended seesaw mechanism. 317 SU(2)l is the usual SM gauge 
symmetry, and Z4 is an additional auxiliary symmetry. For SU(2)l and for A4 the representations 
are given, while for Z4 the charge is indicated. 



Field 


^1,2,3 




I-IR 


TR 


Ni 


N 2 


N a 


Sr 


H 


f 




f" 


5 


e 


X 


SU(2) L 


2 


1 


1 


1 


1 


1 


1 


1 


2 


1 


1 


1 


1 


1 


1 


A 4 


3 


1 


1" 


1' 


1 


1' 


1 


1 


1 


3 


3 


3 


1 


1' 


1 


Z 4 


1 


1 


1 


1 


1 


i 


-1 


i 


1 


1 


i 


-1 


1 


-1 


— i 



Going to neutrinos, the mass matrices can be read off from Eq. (117) as 

( V\ 2/2 
V2 Vi 
. y 2 -2/3 

which yields m s = \p\u as the mass of the keV sterile neutrino. The light neutrino 



m D = — ^- I y 2 y-i I, Mr = &&g(\iv v ,\2V ip , X^vJ), M s = (jm, 0, 0), (119) 



A 



x Note that certain dangerous terms are actually neglected in Ref. 317. It is argued that these 
terms could be forbidden by an additional global 1/(1) symmetry, under which \ an d S are suit- 
ably charged. However, since \ will obtain a VEV later on, this will break the symmetry which 
will unavoidably lead to a massless Goldstone boson. 242, 243 While this problem is not solved in 
Ref. 317, one could at least think about gauging the symmetry or giving a non-perturbative mass 
to the otherwise massless state, making it a pseudo Nambu- Goldstone boson. 
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mass matrix is given by 



m 




(120) 



\ A2 A2A3 A2A3 



which leads to the light neutrino masses 



(mi,m 2 ,TO 3 ) = (0, 



1 2 2 Q 2 2 

' |A 2 |A^ ' |A 3 |A a 



. These 



eigenvalues clearly obey normal ordering. However, the matrix diagonalizing m v 
from Eq. (120) is tri-bimaximal, cf. Eq. (62), which must equal the PMNS matrix 
due to the charged lepton mass matrix being diagonal. Hence, the model presented 
is excluded by the new data on #13, cf. Eqs. (4). It is claimed in Ref. 317 that this 
problem could be cured by next-to-leading order terms, but up to now there exists 
no explicit proof for this claim. 

Nevertheless, other models based on the extended seesaw mechanism could po- 
tentially be found, since there is no principle argument why the extended seesaw 
should be inconsistent with out current data. For example, a very recent proposal 319 
uses anomaly free settings of a U(l)' symmetry instead of A 4 . This has the advan- 
tage to be able to accommodate for the new data, but for the cost of introducing 
some more singlet fields like Sr in order to ensure the cancellation of all anomalies. 
On the other hand, the model has the interesting feature of intrinsically predicting 
a mixture of light and heavy DM. Furthermore, the additional singlets lead to a 
modification of the light neutrino mass matrix such that there exists no massless 
state anymore. Hence, if it turned out that even the lightest neutrino mass is signif- 
icantly different from zero, the {/(l)-extension of the extended seesaw mechanism 
would still have no problems, while its minimal version would be under pressure. 

5.6. More complicated seesaw and suppression scenarios 

Finally, we want to shortly comment on more complicated models which, to some 
extent, rely on a seesaw or similar type of suppression mechanism. These models 
are, from a purely scientific point of view, not inferior to the ones already presented. 
However, for one reason or the other they are slightly more complicated and/or their 
details go beyond the scope of this (hopefully) pedagogical review. Nevertheless, 
these models offer interesting alternatives to motivate a hierarchy in the sterile 
neutrino sector. The reader interested in the finer details is referred to the references 
given in the following. 

5.6.1. Type II seesaw mechanism in a 331-model 

This model has been presented in Ref. 320, based on earlier considerations in 
Ref. 321. It is a slight extension of the ordinary 331-model, 322,323 which extends 
the SM gauge group to a more complicated group SU(3)c x SU(3)l X U(1)n- This 
symmetry is broken at a high scale to the SM gauge group. According to elementary 
group theory, 303 an SU(N) group has N 2 — 1 generators, while a U(l) group has 
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only one. Hence, in addition to the eight gluons from SU(3)c which are unchanged 
in this model, we expect 3 2 — 1 + 1 = 9 gauge bosons in total, four out of which are 
the ordinary electroweak gauge bosons while the other five are new: V ± , U°, U°\ 
and Z'. 

The SM leptons come in SU(3)l triplet and singlet representations, 

1,3,-0 , e«~ (1,1,-1), (121) 



which means that the fields (Nm) c , which we will later on call sterile neutrinos, 
are indeed only sterile under the SM gauge group but not under the full gauge 
symmetry. 

Along with the leptons, we have three triplet scalars and one sextet scalar, 






1,3,--J , (122) 

A° A- *°\ 

1,3, — , S=-= A" A"" $- ~ 1,6, — 
1 1 $o a «) \ 3 

where the $ 8 contain SM-like Higgs components, but they would give a Dirac mass 
to neutrinos. Other SU(2)l doublets, not necessarily listed here, are still available 
to give masses to the other SM fcrmions. 

Since the authors of Rcf. 320 are after a LH seesaw type II mass, cf. Eq. (18), they 
forbid a VEV for $° by imposing suitable discrete symmetries. They furthermore 
tune their scalar potential such that the components (x' , p°, rp , A , a ) obtain 
VEVs (v x ' , v p , v v , «a, v a ). The key point is that some of the scalar components 
do carry lepton number, L{rf ,a a ,p' + ) = —2 and L(x°, X~ ■> A , A - , A ) = +2. 
Hence, the VEV ua (v a ) unavoidably breaks lepton number and leads to a mass 
term for LH (RH) neutrinos. The authors exploit this fact by adding a term to the 
scalar potential which explicitly breaks lepton number, 

Vlnv = -MirfS^r] - Af 2 x T S' t X + (other terms), (123) 

and by furthermore assuming that the corresponding mass scale M sa is larger 
than all VEVs, v x ',p.»7.A.er "C M. This might sound a little surprising since some of 
the VEVs do break lepton number, but in reality one cannot even speak of lepton 
number anymore after it is broken by the explicit violation terms in Eq. (123). 

Moreover, a strong hierarchy between certain VEVs is assumed, Va,o- *C Vx' ,p,n- 
This is only an assumption, but it could be motivated by observing that (A, a) are 
components of the sextet scalar field, while (x'iPiV) are contained in triplets under 
SU(3)l, cf. Eq. (122). This is enforced by the minimization conditions that yield, 
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under the above assumptions, 

v 2 v 2 

m L = Y— , M R = Y^-, (124) 

M M y ' 

which motivates rriz <C Mr. Note that the Yukawa matrices Y are the same for the 
LH and RH fields. 

This is the seesaw-like mechanism had mentioned before: both mass matrices 
in Eq. (124) are suppressed, and in particular we have a motivation for relatively 
small Mr eigenvalues which could well be at the keV scale. However, a problem of 
this model is that it does not explain any hierarchy in the sterile neutrino sector. 
If the goal of model building is to produce settings which resemble the j/MSM, 23 
then this model unfortunately fails to do that job unless even more symmetries are 
imposed. 

Finally note that a variant of this model was presented in Ref. 324, where two 
new SU(3) symmetries, SU(3)l and SU(3)r, appear. In this "3331"-model, each 
of the two SU(3) yields a lepton triplet, which leads to a left-right symmetric 
version of Eq. (121). Without going through the details of this extended model, we 
nevertheless want to mention that this model yields three keV neutrinos Nl which 
have the remarkable property that they do not mix with the ordinary neutrinos, 
thereby invalidating the X-ray constraint, cf. Sec. 4.1. 



5.6.2. Composite Dirac neutrinos 

The other seesaw-type mechanism we want to briefly discuss arises in the context 
of composite neutrinos. Such settings are based on so-called preons, 325 which are 
assumed to be fundamental fcrmions q charged under some new hidden gauge group 
which is often called v-color. Prcons condensate in a QCD-likc manner to bound 
states of a generic mass scale A. The key point is to assume the existence of an 
even more fundamental scale M 3> A, below which everything we talk about is an 
effective field theory. The idea of such a scenario to explain the possible existence 
of keV sterile neutrinos was put forward in Refs. 326, 327. 

All Yukawa couplings of preons to SM- fields are effective vertices, e.g., 

^Yukawa = - M 3(^-l)/2 ^q 71 , (125) 

where A = 0(1) and n is the number of preons q. Certain fields could also be 
charged under some hidden flavor symmetry, which could give a certain structure 
to terms like the one in Eq. (125). Note that, since the n preons obviously have to 
combine to a fermion, and we also need n > 1 to have a bound state in the first 
place, the minimum value possible is n > 3. If the preons condense at some scale 
A to a spin-i bound state, q n ./VrA 3 (™ -1 )/ 2 , then the corresponding Dirac mass 
terms are suppressed by powers of the small quantity e = -jjy , 

^Yukawa -> -XS^^HNr + fl.C, (126) 
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where 3 ("~ 1 ) > 3. Hence, there will be at least an e 3 mass suppression from such 
terms, compared to their natural scale v from the Higgs VEV. However, if also left- 
handed (supposed to be heavy) neutrinos Nl exist, which are also singlets under 
SU(2) - i.e., the fermions N are vector-like - then also other terms of the form 

C LH = -AnNENr + h.c. (127) 

are possible. Here, Ajy = O(A), which would lead to heavier masses. 326 ' 327 As an 
intermediate type of mass term, however, the existence of Nl would also lead to 
masses of the form 

Clw = -A'^'-^NENr + h.c, (128) 

where A^y = 0(A). This would motivate different types of models for the sterile 
neutrinos N, and in particular it could introduce a hierarchy such that some of the 
sterile states have masses of 0(keV). 

The remarkable unique point about this model is that the sterile neutrinos are 
Dirac particles, contrary to all other models presented here: due to the existence 
of the SU(2) singlet fields Nl, one can write down mass terms as in Eqs. (127) 
and (128) for sterile neutrinos, but they will nevertheless be sterile under any SM 
interaction. Furthermore, the mass term in Eq. (126) simultaneously offers an ex- 
planation for the active neutrino masses to be tiny, and the hidden flavor symmetry 
could even be used to obtain a suitable leptonic mixing structure. 

While this model has several nice features, one could object that the new scales 
introduced are completely arbitrary and not even fundamental, since they essen- 
tially rely on some non-perturbativc QCD-likc mechanism. Furthermore, the whole 
mechanism is relatively complicated, which might be another negative point. The 
lesson to learn is that, as always in model building, whenever we want to explain 
something we have to pay a price. 

5.7. Scenarios and models based on other ideas 

Apart from the models discussed in Sees. 5.2 to 5.6, there are several other settings 
which can accommodate for keV sterile neutrinos. However, as we had argued at 
the beginning of Sec. 5, the mere fact that a certain setting can involve keV sterile 
neutrinos, i.e., it does not disagree with their existence, does not a priori yield any 
explanation for them. Hence, in the terminology used in this review, such settings 
would be called scenarios rather than models. 

Let us nevertheless briefly mention a couple of examples from the literature 
which can accommodate for keV sterile neutrinos: 

• keV sterile neutrinos in the scotogenic model : 32S ~ 330 

The scotogenic model is a very simple extension of the SM by three RH 
neutrinos iV,-, one additional Higgs doublet rj, and a Z 2 symmetry under 
which only the new fields are charged non-trivially, i. e. , they receive a factor 
of ( — 1) under this parity transformation. 232 In this setting, there is a very 
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natural explanation for the scale of the light (active) neutrinos, since their 
masses only arise at 1-loop level. This leads to a natural suppression of 
their mass. However, this mechanism to generate active neutrino masses 
intrinsically requires a tree- level Majorana mass for the RH neutrinos A^: 
since the model docs not include any mechanism to spontaneously generate 
any violation of lepton number, which must be present in order for it to 
translate into the light neutrino sector. 

On the other hand, the mass spectrum of the RH neutrinos is not fixed. 
While they could lead to very interesting phenomenology at LHC if their 
masses were of O(TeV) , 328 ' 331-333 it is not excluded that one of the "heavy" 
neutrinos could be as light as O(keV), which is exactly what is considered by 
Refs. 328, 329, 330. While these references discuss very interesting aspects of 
the phenomenology of such settings, they nevertheless lack an explanation 
for the keV scale in the first place. 

• keV sterile neutrinos with Left-Right symmetry: 53 ' 

A very interesting mechanism in order to correct a too large abundance 
is the release of a sufficient amount of entropy by late decay of a non- 
rclativistic species that temporarily dominates the energy density of the 
Universe. 94 This mechanism has also been discussed in the context of 
the ^MSM, where it is able to open up new windows in the parameter 
space. ' 

If a keV "sterile" neutrino is only sterile with respect to SM interactions 
but takes part in gauge interactions beyond the SM, it will naturally be 
produced thermally. 61 However, just as if active neutrinos had a mass of 
O(keV), such a production will overdose the Universe, i.e., the amount of 
DM produced will be too big. This can however be corrected by a sufficient 
entropy production from the decays of the two heavier sterile neutrinos, as 
discussed in Sec. 2.2.4. 

As before, even though this setting can contain keV sterile neutrinos, it 
does not by itself lead to an explanation of their mass. Nevertheless it is 
interesting, in particular since it is not consistent with every way to explain 
the keV scale, cf. Sec. 5.2.1. 

There are also further examples in the literature which would be called models 
according to the terminology used here, but which for one reason or the other do 
not comprise a setting that is suitable for an illustrative example. These include: 

• U(l) symmetries broken close to the Planck scale: 268 

In Rcf. 268, certain structures that can also lead to keV sterile neutrinos 
are discussed. Concretely, a model is presented which is based on a t/(l) x 
U (1)' x Z4 x Z2 symmetry. It connects both the production of the keV sterile 
neutrinos and the generation of their mass to a relatively light inflaton field 
cj). In particular, the continuous symmetries are broken at a scale so high 
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that it could even be close to the Planck scale. 

However, one might object that, in fact, the symmetry used in Ref. 268 
is actually not too much different from a generic Froggatt-Nielsen U(l). 
Furthermore, the breaking of the global continuous symmetry leads to the 
existence of a massless state x, which is on the other hand argued to con- 
tribute to Dark Radiation 208 ' 334,335 (see Ref. 54 for an extensive collection 
of oservations) - this could even be desirable from a cosmological point of 
view. On the other hand, certain potentially dangerous terms are neglected 
in the scalar potential, and the above mentioned Z 2 symmetry is completely 
artificial as it has no effect other than stabilizing the lightest sterile neu- 
trino Ni and by this avoiding the X-ray bound, cf. Sec. 4.1, in order to 
make the model consistent with current (non-) observations. 

• light sterile neutrinos from global symmetries: 336 

Another approach is to generate the sterile neutrino mass only by an 
effective operator, which itself arises from an unknown (or partially un- 
known) high energy sector. This possibility has been extensively discussed 
in Ref. 336, where exactly one sterile neutrino called vg is introduced, which 
receives mass from a Weinberg-like operator. The authors mention the pos- 
sibility of an operator \{vs) ClJ sHiH2 with two Higgs doublets #1,2, but 
then discard this possibility because they are after 0(l)-active sterile mix- 
ing with a sterile neutrino mass of about 1 eV, in order to explain the by 
then famous LSND anomaly. 337 Instead, they use an operator j^(vs) c vsSS 
with a two singlet scalars S. While a bare (and hence potentially large) mass 

term msiys) Cv S is forbidden by a global symmetry, the VEV (S) = vs rein- 

2 

troduces the possibility of a higher order mass ^f- for the sterile neutrino. 
In particular, the authors discuss how to choose the required global sym- 
metry appropriately without being in danger to end up with anomalies, cf. 
Ref. 336 for details. 

While the general approach of Ref. 336 is certainly valuable, the reference 
lacks the details of the high energy sector, even though the results derived 
are not very sensitive to that. However, a more serious drawback is that 
there is no actual spectrum present in the sterile sector, so that the resulting 
neutrino mass is only a function of two unknown scales, vs and A. Never- 
theless, the authors approach certainly fulfills the criterion to completely 
disentangle the mass generation of the sterile neutrino from the electroweak 
scale, which allows for this freedom. 

• keV sterile neutrinos from Dark GUTs: 338 

The possibility of a seesaw mechanism could not only be used in an active 
sector, but also in a dark sector of a theory which hardly talks to the SM 
part. Such an approach has been taken in Ref. 338, where different groups 
G' are studied to be used in a setting similar to a Grand Unified Theory 
(GUT), which however does not include the SM itself as a gauge group. The 
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field which plays the role of the sterile neutrino is charged non-trivially 
under this new symmetry G", and in this sector it receives a suppressed 
mass in pretty much the same way as an active neutrino receives its mass 
in the ordinary seesaw type I mechanism, cf. Appendix A. This field is 
nevertheless a singlet under the SM, and it can hence play the role of a 
sterile neutrino. 

While a couple of choices for the group G' are discussed in Ref. 338, the ap- 
proach is nevertheless somewhat sketchy Furthermore, although a general 
argument for the smallness of the sterile neutrino mass scale is given, there 
is no reason for enforcing a certain structure in the sterile sector, which 
could lead to phcnomcnological problems. 

• loop-masses for keV sterile neutrinos: 339 

A very natural mechanism to explain hierarchies in the sterile neutrino 
sector has been presented by Ma in Ref. 339, by using loop suppressions of 
the masses of certain sterile neutrinos. In the model suggested, SM singlet 
fcrmions 51,2.3 arc introduced in addition to the "ordinary" RH neutrinos 
A^i.2,3, along with several additional scalar fields. The arrangement of the 
fields and charges is such that a beautiful thing happens in the sterile sector: 
the field Si is the only additional singlet which has a tree- level mass Msi- S2 
in turn, only receives a mass at 1-loop order by a diagram which involves 
the massive fermion Si in the loop, leading to a suppressed mass M$2- 
Finally, the field S3 also receives a mass by a 1-loop diagram that contains 
the field S2 as internal fermion, which by itself only receives a 1-loop mass. 
This means that the corresponding mass M$3 is in fact suppressed by a 2- 
loop mechanism. Such a setting naturally yields an explanation for M53 <C 
Ms2 *C Msi, and by this it is a good motivation for M$y, = O(keV). 
However, what had been presented above is the beauty, and the beast comes 
along in the actual charge assignments: even though the above mechanism 
is plain and simple, the charge assignments necessary to generate it are, 
unfortunately, very complicated. Not only does the model require many 
more or less unmotivated scalar fields, it also requires an additional U(l) 
symmetry with either a very large charge of one scalar or with relatively 
tricky charge assignments for the singlet fields 5*1,2,3- While this model could 
certainly yield interesting phenomenology, its setting to some extend looks 
artificial and engineered, which is the price one has to pay for the elegant 
mechanism illustrated in the previous paragraph. 

• light sterile neutrinos in minimal radiative inverse seesaw: 340 

The setting of Ref. 340 is based on the so-called inverse seesaw 230,341 
and it is, in effect, a combination of the seesaw mechanism with a loop- 
suppression. The idea is to extend the seesaw type I setting by further 
LH singlet fermions Sl, which themselves obtain a very large Majorana- 
type mass /is and which also form a Dirac mass term with the ordinary 
right-handed neutrinos Nr. Integrating out the heaviest singlets Sl then 
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leads, at tree-level, to a mass matrix that looks very similar to an ordinary 
type II mass matrix, and in particular it has a seesaw structure in all en- 
tries. However, it turns out that the rank of this matrix is nevertheless too 
small to yield light sterile neutrino masses (in other words, a cancellation 
forces them to be zero). On the other hand, at 1-loop level this mass ma- 
trix receives a correction which induces small sterile neutrino masses that 
are naturally loop-suppressed in addition to the seesaw structure that was 
already present. One could say that the loop-correction destroys the exact 
cancellation just mentioned. 

This mechanism is very strong, so that sterile neutrino masses of O(keV), 
or even lighter, are no problem to achieve. Furthermore, this mechanism is 
also quite minimal, in the sense that the seesaw framework is not extended 
by very much, while the suppressions follow naturally from the general 
structure. However, one might ask why the Majorana neutrino mass for the 
left-handed singlets is so much bigger than the (equally unconstrained) one 
for the right-handed neutrinos. In addition, a further potential issue is that 
this mechanism again leads to a very strong suppression of certain sterile 
neutrino masses, but there is no argument for a particular mass pattern. 
• anomalous sterile neutrino masses from gravitational torsions: 342 

A more exotic way of generating a light sterile neutrino mass is presented 
in Ref. 342. The mechanism relies on so-called quantum torsions, which 
may arise in string theory or in more general quantum gravity settings. 
The general idea is to introduce two axion-like pseudoscalar fields a and 
6, one providing the chirality change in the coupling of two RH neutrino 
fields Nr and one being non-propagating and providing the coupling to the 
quantum torsion. Such a scenario allows to draw a two-loop diagram giving 
mass to Nr, which involves the pseudoscalar a, linearized graviton fields, 
and a coupling to the torsion. 

The setting yields two possible suppressions: with only one axion a, the 
diagram is strongly UV-divergent and the mass Mr is proportional to the 
sixth power of the cut-off scale A. Hence, for a low string mass/cut-off 
close to the GUT scale, Mr can be of O(10 keV) for a Yukawa coupling 
y a = 0(1O~ 3 ) in the term proportional to o,(Nr) c Nr. Alternatively, if there 
are at least three axions a n (n > 3), the mass Mr can be made finite (at 
two- loop level), with its size being controlled by the ratio of axion mass- 
splittings to the axion mass scale, ( -jjt ) ■ 



Although this mechanism yields a relatively general suppression for geome- 
tries with torsion in a quantum gravity setting, one could criticize that it 
relies on a relatively complicated framework, and hardly any scales and 
couplings in the game are fixed or at least tied to some other scale. Such 
ties could be naturally given in a more concrete model which, however, does 
not exist up to know - or at least the author is not aware of it to the best 
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of his knowledge. 

Let us end this chapter by commenting on a generalization of keV sterile neutri- 
nos, and other quasi-sterile fermions. Extremely weakly interacting particles with 
a mass of around a few keV acting as (typically warm) Dark Matter can also be 
found in frameworks other than keV sterile neutrinos, for example axinos, 343 grav- 
itinos, 343-345 KK-gravitons, 343 majorons, 346 ' 347 mixed axion-axino DM, 348 moduli- 

nog 349-351 Qr ginglinos. 352 

This observation has led the authors of Ref. 97 to a generalization of the 
fermionic WDM candidates to a new class of DM called ke Vins, which stands for 
keV inert particles. Just as WIMPs (Weakly Interacting Massive Particles), keVins 
comprise a relatively general class of DM particles, of which keV sterile neutrinos 
are only one special case. However, also for keVins there is in general no reason 
to have a mass of a few keV, just as there is no reason for a general WIMP to be 
weakly interacting. Such a reason has to be imposed in a concrete model. 

Note that in Ref. 97, the mechanism to produce the correct DM abundance is, 
just as in the LR-symmetric versions of keV sterile neutrinos, 53 ' 91 thermal overpro- 
duction with subsequent dilution by entropy production. 94 For keVins, the particle 
that decays and produces the required entropy is a heavier version of the DM particle 
with a mass of a few GeV: the GeVin. However, this method is only one possibility 
to achieve the correct DM abundance, and even the production mechanism men- 
tioned here has not yet been studied in very great detail. In general, the class of 
keVins offers more possibilities for relatively general studies, which could then be 
used to discover new successful regions in the parameter space which might not be 
visible in all concrete models. 

In this central chapter of the review, we have seen that the key point in keV 
sterile neutrino model building is to find a motivation for the keV scale in the 
first place. This by itself is far from trivial, and most of the known mechanisms 
achieving an explanation have been presented in Sees. 5.2 to 5.6. Interesting further 
possibilities arise when trying to generate keV neutrino masses, e.g. at loop-level, 
but although several such settings are available for light neutrinos, there is still no 
simple model to generate a sterile neutrino mass in that way. Of course, further 
mechanisms and arguments could be found, and by this we could not only establish 
links between two very fundamental sectors in particle physics, neutrinos and Dark 
Matter, but we could also gain an improved understanding of our model building 
techniques by investigating in how far such an extremely hierarchical mass spectrum 
as needed for keV sterile neutrinos could be realized. 

6. Conclusions 

We have presented a review of the current ideas present in the literature which can 
explain or at least motivate the existence of a keV sterile neutrino. Such a particle 
would be well suited to act as (typically warm) Dark Matter and they can yield 
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very interesting connections to the low energy neutrino sector. While keV sterile 
neutrinos are very well investigated from the astrophysics side by studies of their 
production mechanisms and their behavior in cosmology, the model building side has 
only become a rising field within the last few years. Nevertheless, a review at exactly 
this stage is useful in order to motivate the scientists from the different disciplines to 
collaborate on this exiting topic. It has become clear for model builders that there 
are many interesting relations of keV sterile neutrinos to the light neutrino sector, 
which could lead to both, new types of mass-creating mechanisms and testability 
by the light neutrino sector. These interrelations are under investigation by several 
groups by now, but nevertheless the corresponding scientific works have sometimes 
been hard to understand for people from neighboring fields. This gap has hopefully 
been bridged at least to some extent by the present review. 

Ideally, we should soon enter an era where model studies are combined with 
detailed studies of the cosmological implications of sterile neutrinos. Many exciting 
new experiments and observations are on the way, and they can be expected to 
have a considerable impact on the field within the coming years. We should take 
the opportunity we have now to be prepared in order to make optimal use of the 
data to be gained. New discoveries are the advances of science, and the growing 
field of keV sterile neutrino Dark Matter will hopefully benefit from several such 
discoveries in the near future. 
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Appendix A. The seesaw mechanism in detail 

In this appendix, we perform the seesaw calculation in great detail. Although known 
in principle, in many references (sometimes even in the supposed-to-be-pedagogical 
ones) this calculation is not presented in sufficient detail, and certain logical steps 
are skipped over. Due to its importance, we will go through it without skipping 
anything. We will furthermore explain two related technicalities which are of interest 
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here, namely the so-called Casas-Ibarra parametrization and the notion of form 
dominance. 

Appendix A.l. The seesaw mass matrix 

The Dirac mass term for neutrinos is given by 

£d = -Wm D N R + h.c, (A.l) 
while the LH and RH Majorana mass terms are given by 

1 , 1- 



C M = -^Em L {iy L ) c - -{N R ) C M R N R + h.c. (A.2) 

The combination of the two terms, C v = Cu + Cm , leads to the well-known gen- 
eral seesaw type I 220 ~ 224 and type II situations, 226,227 while the limiting case where 
rriL = is just the type I seesaw and the limiting case where tud = Mr = 0, where 
rriL arises due to a Higgs triplet field, is just type II. Note also the structural dif- 
ference between the LH and RH Majorana mass terms in Eq. (A.2): one could have 
defined them in an analogous matter, which would simply amount to a redefinition 
of either mi or M R , but the choice made here will turn out convenient later on. 

Note that Majorana mass matrices are symmetric, = and Mj£ = M R : 
for any Majorana spinor we have $ c = e"^^ , with a Majorana phase </>, from 
which it follows that 

=i 

= e-^ T p CMj a CWf = e^"-W p Mjj>% = W p M$Ji a , (A.3) 
where the first minus sign arises from the anti-commutation of fermion fields. We 

— T 

have also used the ordinary charge conjugation formula w c = Cw , which holds for 
any spinor ^> , Majorana or not, with the charge conjugation matrix C = ij 2 j° . We 
have further exploited that 

f c = C* T = C(* t 7 °) T = C*7° =S> 7 °* c = -C \F = * T C. (A.4) 

- 7 °C Ct 

Renaming the indices on the right-hand side of Eq. (A.3) finally leads to 

^MapVp^MZpVp, (A.5) 

from which we can immediately read off that M = M T . 

The final trick to arrive at the correct form of the seesaw mass matrix is to split 
up the Dirac mass term, Eq. (A.l), into two pieces: 

Wm D N R = ^T7Lm D N R + \ ) (yLm D N R ) T . (A.6) 
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Making use of the anti-commutation of fcrmion fields as well as of Eq. (A. 4) one 
more time, we can rewrite the transposed term: 

{Wm D N R ) T = -(N R ) T m T D (~C 2 )W T = (N R ) T C m T D CW^. (A.7) 
Putting all together, we finally obtain: 

C v = C D + C M = --Wm D N R ~ ^{N R ) c m T D {v L ) c - ^T7j;m L {u L ) c - ^(N R ) C M R N R + h.c. 
m L m D \ ( (u L ) 



=- 5 (-.<w^m;j^J + ^ (a - 

This equation involves the full (typically 6x6) neutrino mass matrix. 
Appendix A. 2. The seesaw diagonalization 

The next step is to (approximately) diagonalize the full neutrino mass matrix: 

M v = ( mL T 2 D ) , (A.9) 
\rrip M R J 

where one assumes that somehow tul <C mo <C M R . Of course, this is an intrinsi- 
cally ill-defined condition when dealing with matrices, in particular since even the 
"large" matrix M R could have texture zeros. Still, we will follow the practical as- 
sumption that the hierarchy between the mass matrices is justified, but we should 
keep in mind that in some cases it may be useful to numerically check if the low 
energy mass matrix actually does conserve the structure of the full mass matrix. 

Similar diagonalization procedures as the one here are outlined in, e.g., Refs. 273, 
353, although we slightly depart from both references. The first step is to approxi- 
mately block-diagonalize the matrix by a nearly unitary matrix U given by 

u ={ -h i-p b )> (AJ0) 

where the parameter b, to be defined later, is assumed to be of OIttidM^ 1 ), and it 
can hence be regarded as small. In particular, since is the smallest mass scale 
in the game and since we expect 0(t7il) = Q(m 2 D M R ~ ) for type II seesaw, we will 
only need to keep mj, with a coefficient of 0(1), mo up to a coefficient of 0(b), 
and Mr up to a coefficient of 0(b 2 ). Then, we can immediately see that 

W-^-( 1 + * " W 1 + ^ lsl4 )-l + ^)-l. (A.U) 

The next step is to calculate D = U T M v l] , keeping only relevant (i.e., large enough) 
terms. For example, we have 

rriL — ^b*b'mL — b*rnjj mo — \b*b T mu — b*Mf 
i T D -±b T b*m T D ' 

/m/, — 6*77Tq mr) — b*M R 

\ m T D b T m D + M R - \b T b*M R 



b T m L + m T D - \b T b*m T D b T m D + Mr - \b T b*M R 



(A.12) 
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All that remains is to calculate the product, 

~ (m L -b*m T D 77i d — b* Mr \(l-\bW b \ , 

V rn T D b T m D + M R ~\b T b*M R )\ -ftt l-\tfb)' [AA6) 

The off-diagonal elements of this expression should vanish, which determines the 
matrix b unambiguously. For the 12-element, we have 

D 12 ~ m L b-b*mnb+m D -b*M R ~ m D -b*M R = &* = moM^ 1 , (A.14) 

which leads to b = m* D M~^ , b T = M^ 1 m^, and = b T * = M^m^. This also 
implies that D 21 — 0. Similarly, we can calculate the diagonal blocks: 

D\i ~ rtiL — b*7nj) — mptf + b*M R ¥ = rriL — mjjMZ mi 



i D — iiojju t u i\± R u — iiol — ub]ji\i R nor), 

D 22 ~M R +- 



mlmoMj^ 1 * + M R l *m^ D m D ~M fi . (A.15) 
By this we have recovered the general seesaw formula, 

D ee U T M U U * ( mL - m ° M R m n ) = ( m » ) . (A .16) 
V M RJ \ M RJ 

In the basis where the charged lepton masses are diagonal, the remaining 3x3 
blocks can be diagonalized by C/p MNS m v Z7pMNS = diag(mi, m 2 , ms) = D v and 
V R M R V R = diag(Mi, M 2 ,M 3 ) = D N , which finally leads to 

D = U T MJI= ( Dv ° j , where (A.17) 



. D N/ 

jj (l - \m* D M^* M^ml m* D M^* \ (U PM ns 

~ \ -M^ml 1 - ^M^mlm^Mn 1 * J \ V R 

which is a central result of this appendix. 



Appendix A. 3. The Casas-Ibarra parametrization 

In the context of the seesaw type I mechanism, a smart parametrization has been 
proposed by Casas and Ibarra in Ref. 354. The trick of this parametrization is to 
make maximal use of the observable parameters in a seesaw type I setting, while all 
unknowns are cast in an orthogonal matrix R. We will shortly review the derivation 
of the parametrization in the conventions used in this review. 

The starting point is the type I seesaw formula for light neutrinos, cf. Eq. (A. 16), 

m v = — to,dM^ 1 to5, (A. 18) 

in a basis where the charged lepton and RH neutrino mass matrices are diagonal. 
Then, exactly as in Eq. (A.17), m v can be diagonalized by the PMNS matrix U = 
f^PMNS, D„ = diag(mi, m 2 , m.3) = U T m v U 1 where U has exactly the form given in 
Eq. (3). With the following definitions, 

ee diag( v /roT, ^/mj, y/frij) and D^L = diag(l/ y/~M\, 1/ \fM 2 , 1/ \fM^), 

(A.19) 
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one can rewrite the diagonalization of Eq. (A. 18), 

D^D^ = D V = U T m u U = -U T m D M^m T D U = U T im D D^LD^LimlU. 



(A.20) 

Multiplying from left and right by D~^. one obtains 

D^LimlUD^L] T [D-^imlUD-^] , (A.21) 

which suggests the definition of the orthogonal complex matrix 

R = iD^LmlUD^. (A.22) 

Note that the above parametrization can also be extended to seesaw type II 
scenarios: 355 the trick is simply to make the replacement 

m„ — > X v = m v — rriL (A. 23) 

in Eq. (A. 18), where tul is the LH neutrino mass matrix, cf. Eq. (A. 8). Of course, 
this also means that in the subsequent equations one has to replace 

D u ^D x =di&g{X 1 ,X 2 ,X 3 ) and U^U X , (A.24) 

where Xi are the eigenvalues of the matrix X v and Ux is the matrix that diagonal- 
izes X v according to Dx = diag(Ai, X2, A3) = V^XyU x- While by this one loses 
the direct connection to the PMNS matrix, it is nevertheless useful since one can 
still exploit the general parametrization of the complex orthogonal matrix R, 

R = i?i 2 i?i 3 i? 23 , (A.25) 
where R^ are the rotations by complex angles cjy = pij + iutj in the ij-plane, 

Su 13 

R12 = I s ui2 c wi , I , R13 =1 10 




where c Uij = cosojy and s u 







Rl3 = ^ 









= sinoj,;j 





0^13, 





f 1 







R23 — 







') 

°(^23 1 











(A.26) 



Appendix A. 4. Form dominance 

The notion of form dominance was defined in Rcf. 260. In a seesaw type I situation, 
we first denote the columns of the Dirac mass matrix by 3- vectors in a basis where 
the right-handed neutrino mass matrix is diagonal/ 

m D = (A, B, C), Mr = diag(A/ A , M B , M c ). (A.27) 
Using this notation, the seesaw formula can be easily shown to have the form, 

T (a®a t b®b t c®c t \ , aoo , 

m v = -mnM^ml = - + + , (A.28) 



y Only here, we change the notation of the right-handed masses to Ma,b,Ci since it will make it 
easier to follow the calculation 
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where the tensor product is defined in the usual way, e.g., 

I MM A,A 2 A 1 A 3 \ /AA 
A® A 7 = A 2 A X A 2 A 2 A 2 A 3 , where 1= \A 2 \. (A.29) 
\A 3 Ai A 3 A 2 A 3 A 3 J \A 3 J 

Let us now write the PMNS matrix U in terms of its column vectors, 

U = {Ui,U 2 ,U 3 ). (A.30) 

Then, in the basis where the charged lepton mass matrix is diagonal, form domi- 
nance is defined as the situation where each column of mr> is proportional to the 
corresponding column of the conjugate PMNS matrix U* , z 

form dominance (A, B, C) = (aU{,bU 2 ,cU 3 ), (A.31) 

where a, b, c are arbitrary complex numbers. Note that, once a proportionality of 
the columns is given, we can always write this as in Eq. (A.31), by relabeling the 
columns A, B, and C . Explicitly, one can write: 

(U el \ (U e2 \ (U e3 \ 

A = a [U^ , B = b C^a , C = c [U^ . (A.32) 
\UrJ \U T2 J \U t3 J 

Since the columns of a unitary matrix are orthonormal, Ujllj = 5ij and tjftj* = <5, j, 
one can easily show that: 

/l 0\ /0 0\ /0 0\ 

U T A®I r U = a 2 , U T B®B T U = b 2 1 , U T C®C T U = c 2 . 

\0 0/ \0 0/ \ooi/ 

(A.33) 

This allows to express the light neutrino masses in terms of the heavy neutrino 
masses and only three free parameters, 

a 2 b 2 

Ma ' _ Mb ' _ McJ 
(A.34) 

Finally, we note that form dominance is equivalent to setting R = 1 in the Casas- 
Ibarra paramctrization, cf. Appendix A. 3. 356 Indeed, using Eq. (A. 22), one obtains 

m D = U*D^{-iR T )D^ = U*D^ = -i^rmM^, ^m 2 M 2 U* 2 , ^/m 3 M 3 U;), 

(A.35) 

which just coincides with Eq. (A.31). 



rrT rr rrT A® & B ® B C®C T \ TT / 

U 1 m v U = -U 1 — — 1 — H — U = diag 

I M A M B M c j V 



z Note that in Ref. 260 U was taken to be real, which is why the complex conjugation does not 
appear in the corresponding formula. 
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